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ABSTRACT 


The  thickness  of  the  Athabaska  glacier  has  been 
obtained  along  eight  transverse  profiles  by  an  investiga¬ 
tion  of  gravity  anomalies.  Three  dimensional  computations 
with  a  low  speed  digital  computer  were  made  in  this 

5' 

study  to  acquire  more  precise  results  than  previously 
obtained.  The  depth  of  the  glacier  varied  from  1070  feet 
on  a  line  below  the  lower  ice  fall  to  160  feet  near  the 
terminus.  The  accuracy  of  the  results  are  discussed  and 
compared  with  independent  data  from  boreholes  and  a 
seismic  program.  The  cross  section  of  the  valley  of 
the  glacier  was  found  to  approximate  a  parabola  on 
several  lines.  From  a  knowledge  of  the  depth,  shape  and 
surface  slope  of  the  glacier,  the  average  shear  stress 
exerted  by  the  bed  on  the  ice  was  found  to  be  1.0  bars. 

Gravity  studies  along  the  Rocky  Mountain  Trench 
south  of  Cranbrook,  British  Columbia,  revealed  the 
presence  of  three  deep  basins  filled  with  unconsolidated 
material.  The  basins  range  in  depth  from  1400  to  3500 
feet  (-17$  to  +25 ft).  Combined  gravity  and  geologic 
data  indicate  that  the  Trench  south  of  Fort  Steele  is 


due  fundamentally  to  block  faulting  with  the  amplitude 
of  faulting  increasing  towards  the  International  Boundary. 
A  linear  gravity  minimum  is  associated  with  the  Moyie 
Lenia  fault  southwest  of  Cranbrook.  This  trend  is 
really  composed  of  two  high  anomalies,  one  on  each  side 
of  the  fault,  which  are  caused  by  dense  sills  of  basaltic 
material  intruding  the  Aldridge  formation.  The  upthrown 
block,  on  the  northwest  side,  brings  the  Aldridge 
formation  and  sills  to  the  surface.  From  geologic 
evidence  the  down  thrown  side  has  a  vertical  displace¬ 
ment  of  about  15,000  feet.  To  satisfy  the  gravity  data 
the  total  thickness  of  the  sills  and  intervening  quart¬ 
zites  should  increase  from  2000  to  5000  feet  as  one 
passes  over  the  fault  towards  the  southeast.  Gravity 
calculations  were  made  treating  the  break  as  a  steeply 
dipping  overthrust  and  as  a  normal  fault.  Of  the  two 
interpretations  presented  the  normal  fault  solution  fits 
the  observed  gravity  data  better.  Important  faults  such 
as  the  Dibble  Creek  and  Boulder  Creek  faults,  which 
appear  in  the  Rocky  Mountains,  extend  at  least  four 
miles  Into  the  trench. 
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1.  INTRODUCTION 


In  the  summer  of  1959  the  University  of  Alberta  and 
the  University  of  British  Columbia  undertook  an  expedition 
to  the  Athabaska  Glacier  to  study  particular  problems  in 
the  geophysics  of  glaciers.  The  program  was  sponsored 
by  the  National  Research  Council  of  Canada  and  was  under 
the  direct  supervision  of  Professor  J.  A.  Jacobs  and 
Professor  G.  D.  Garland.  The  major  factor  determining 
the  selection  of  the  glacier  to  be  studied  was  ease  of 
access.  The  Athabaska  Glacier  was  an  obvious  choice  as 
it  is  a  smooth  valley  glacier  with  a  well  defined  terminus 
which  is  only  one -half  mile  from  a  major  paved  route. 

It  Is  also  65  miles  from  the  town  of  Jasper  on  the  main 
railroad  line  between  Edmonton  and  Vancouver.  Heavy- 
equipment  was  transported  two  and  one -half  miles  from 
the  toe  of  the  glacier  to  the  7300  foot  level  by  Snow¬ 
mobiles  operated  by  Mr.  Wm.  Ruddy. 

The  major  part  of  the  program  was  designed  to  deter¬ 
mine  the  mechanics  of  glacier  flow.  In  this  connection 
80  stakes  were  set  out  for  a  3tudy  of  surface  movement 
and  ablation  by  Stan  Paterson  of  the  University  of 
British  Columbia.  The  variation  of  velocity  with  depth 
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was  to  be  determined  by  long  tern  study  of  inclinometer 
measurements  in  cased  boreholes  following  a  method 
pioneered  by  Gerrard,  Perutz  and  Roch  (1952).  John 
Stacey  was  in  charge  of  designing  the  hot points,  the 
drilling  procedure,  and  general  camp  routine. 

The  thickness  of  the  ice  along  a  longitudinal  tra¬ 
verse  was  determined  by  reflection  seismograph  methods 
with  a  set  of  Houston  Technical  Laboratories  High  Resolu¬ 
tion  instruments.  One  borehole  penetrated  the  ice 
column  to  bedrock  and  provided  valuable  velocity  control 
for  this  seismic  study.  Professor  J.  C.  Savage  was  in 
charge  of  the  seismic  operations.  The  thickness  of  the 
ice  and  the  shape  of  the  bedrock  along  transverse  profiles 
was  to  be  obtained  by  a  gravity  investigation  using  the 
Worden  Gravity  Meter.  The  gravity  study  was  planned  to 
provide  a  detailed  description  of  the  thickness  of  the 
glacier  so  that  the  movement  studies  could  be  properly 
evaluated  and  the  average  shear  stress  on  the  bedrock 
calculated.  It  would  then  form  a  valuable  adjunct  to 
the  seismic  program  which  was  much  more  expensive  to 
carry  out  and,  therefore,  limited  in  scope.  In  addition, 
the  gravity  data  checked  areas  with  poor  seismic  reflec¬ 
tions  and  depth  points  determined  from  events  suspected 
of  being  multiple  in  origin. 


. 
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In  addition  to  the  above  studies.  Professor  G.  D. 
Garland  made  temperature  measurements  In  uncased  boreholes 
to  determine  the  thermal  regime  of  the  glacier.  Inde¬ 
pendent  investigations  were  carried  out  by  other  parties 
on  the  3ame  glacier  during  the  1959  season.  These 
included  a  detailed  aerial  photographic  and  topographic 
survey  by  the  Water  Resources  Branch  of  the  Department 
of  Northern  Affairs  and  National  Resources.  This  survey 
established  a  series  of  .valuable  benchmarks  around  the 
perimeter  of  the  Athabaska  Glacier.  Various  experiments 
to  determine  the  thickness  and  electrical  properties  of 
an  ice  sheet  by  electrical  and  electromagnetic  methods 
were  conducted  by  A.  Don  Watt  and  Gene  Maxwell  of  the 
United  States  Bureau  of  Standards  and  by  Henri  Vetter 
from  the  University  of  Alberta. 

The  gravity  investigation  of  the  glacier  forms  the 
first  topic  of  this  thesis.  The  second  part  of  the 
thesis  is  a  gravity  investigation  of  the  Rocky  Mountain 
Trench  between  Cranbrook,  British  Columbia,  and  the 
International  Boundary.  The  two  topics  are  related 
because  the  major  anomalies  in  both  areas  are  prominent 
negative  minima.  At  the  Athabaska  Glacier,  the  anomaly 
is  due  to  low  density  ice  filling  a  mountain  valley 
while  in  the  Trench,  the  anomalies  are  due  to  depressions 


. 
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in  the  bedrock  filled  with  low  density  unconsolidated 
material.  In  both  cases  one  face  of  the  body  producing 
the  minimum,  that  is  the  surface,  is  completely  determined. 
The  interpretation  of  the  gravity  data  is,  therefore, 
similar  in  principle. 
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PART  I 

AN  INTERPRETATION  OF  SOME  GRAVITY  MEASUREMENTS 
ON  THE  ATHABASKA  GLACIER,  ALBERTA 
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2.  PHYSICAL  DESCRIPTION  OP  THE  AREA 

2. 1  The  Columbia  Icefield 

The  Athabaska  Glacier  is  the  second  largest  ice 
stream  issuing  from  the  Columbia  Icefield,  an  ice  cap 
straddling  the  Continental  Divide  which  marks  the  boundary 
between  the  provinces  of  Alberta  and  British  Columbia. 

The  Columbia  Icefield  covers  an  area  of  110  square  miles 
and  is  contained  by  a  series  of  peaks,  nine  of  which  are 
over  11,000  feet  above  sea  level.  These  include 
Mb.  Columbia  (12,294'),  and  the  North  Twin  (12, 085* ), 
respectively  the  second  and  third  highest  peaks  in  the 
Rocky  Mountains.  The  area  of  accumulation, with  an 
average  elevation  of  nearly  10,000  feet  and  a  maximum  of 
11,340  feet  on  the  Snow  Dome, has  an  abundance  of  snow 
precipitation  in  the  winter  and  throughout  the  year. 

The  elevation  is  sufficiently  large  for  the  temperature 
to  be  low  so  that  much  of  the  fallen  snow  is  conserved 
and  transformed  into  ice.  (See  Figure  1  and  Plate  I.) 

The  Columbia  Icefield  and  Athabaska  Glacier  were 
discovered  by  Professor  J.  N.  Collie,  H.  E.  M.  Stutfield 
and  Hermann  Wooley  in  1897*  In  1919  Dr.  A.  0.  Wheeler 
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PLATE  I.  Two  photographs  of  the  Athabaska  Glacier  from  Mount 
Wilcox.  Mount  Athabaska  (11,432')  is  on  the  left  side.  The 
lower  photograph  shows  details  of  the  three  steps  which  lead 
to  the  Columbia  Icefield.  Enow  Dome  (11,340')  is  out  of 
sight  on  the  upper  right  hand  side. 
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and  R.  W.  Cautley  led  parties  which  surveyed  and  photo¬ 
graphed  the  region  of  the  Continental  Divide  for  the 
Canadian  Government.  The  topographic  maps  produced  by 
the  Interprovincial  Boundary  Commission  are  on  a  scale 
of  1:62,500  and  have  a  contour  interval  of  100  feet. 

(See  Figure  1.)  These  maps  have  proven  invaluable  in 
computing  the  topographic  correction  for  this  gravity 
survey . 

2.2  The  Athabaska  Glacier 

The  glacier  Issues  from  the  Columbia  Icefields  at 
an  elevation  of  about  9000  feet  and  flows  4  miles  to  the 
northeast  to  a  well  defined  terminus  at  6320  feet  above 
sea  level.  The  upper  half  of  the  glacier  descends  in 
three  sharp  drops  to  form  two  prominent  steps.  The 
lower  half  is  2-1/4  miles  long  and  descends  smoothly 
from  an  elevation  of  7500  feet  with  a  dip  of  3-1/2°  at 
station  L  10,  decreasing  to  2-1/2°  at  L  22  and  gradually 
increasing  to  10°  at  L  36  and  about  30°  at  the  terminus. 
The  majority  of  the  observations  were  made  on  this  smooth 
lower  half.  The  fresh  winter  snow  retreats  from  the 
lower  part  of  the  glacier  during  the  early  part  of  July 
and  towards  the  end  of  the  month  the  ice  is  exposed  all 
the  way  to  the  second  step.  The  firn  line  is  on  the 
upper  part  of  the  ice  drop  from  the  icefield.  The  glacier 
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is  the  source  of  the  Sunwapta  river  which  joins  the 
Athabaska,  Slave  and  Mackenzie  rivers  to  flow  2200  miles 
into  the  Arctic  Ocean. 

European  records  of  glacier  retreats  and  advances 

are  fairly  complete  as  far  back  as  1850  when  a  pronounced 

advance  was  recorded  in  the  Alps,  Iceland,  and  Greenland. 

A  recession  was  recorded  between  i860  and  1870  and  a 

minor  advance  between  1875  and  I89O.  Since  the  1890' s 

there  has  been  a  slow  general  recession  in  Europe  and 

America.  Regular  observations  on  the  terminus  of  the 

Athabaska  glacier  were  begun  in  1945  by  the  Water 

Resources  Branch,  Department  of  Northern  Affairs  and 

National  Resources.  The  amount  of  recession  between 

1945  and  1956  amounted  to  468  feet,  which  is  an  annual 

recession  of  58.5  feet.  During  the  same  period  a  series 

of  plaques  on  the  snout  of  the  glacier  advanced  at  a 

mean  annual  rate  of  51  feet.  (Collier,  1957)* 

Additional  data  can  be  obtained  from  the  1919 

Boundary  survey  maps  and  photographs  made  by  Cautley  and 

Wheeler.  These  indicate  a  recession  of  2700  +  200  feet 

•% 

in  40  years  or  an  annual  mean  recession  of  68  +  5  feet. 
2.3  Geology  of  the  Area 

The  Athabaska  glacier  has  formed  a  terminal  moraine 
nearly  3/4  mile  long  during  its  retreat.  On  either  side 
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it  is  bounded  by  high  lateral  moraine  made  up  of  loose, 
fairly  fresh  looking  limestone  debris.  To  the  southeast 
the  lateral  moraine  is  70  feet  above  the  ice  surface  and 
extends  from  the  E  line  to  the  terminus.  It  is  not 
known  if  ice  underlies  this  material.  The  wide  northwest 
embankment  is  covered  with  loose  and  unstable  slabs  of 
limestone  rock  and  is  known  to  contain  glacier  Ice  under 
the  material.  The  ice  outcrops  in  several  places, 
notably  between  lines  H,  F  and  G.  The  known  exposures 
are  30  to  70  feet  above  the  level  of  the  central  clean 
ice  and  is  protected  from  ablation  by  the  mantle  of 
rock.  An  alluvial  plain  extends  along  the  Sunwapta 
valley  for  several  miles. 

No  detailed  geologic  work  has  been  published  for 
this  area.  J.  A.  Allan  (1938)  made  a  regional  survey 
of  the  immediate  area  and  R.  D.  Hughes  (1953)  has  made 
a  detailed  study  of  the  Sunwapta  Peak  area  some  10  miles 
to  the  north.  J.  L.  Severson  (1950)  studied  the  Devonian 
stratigraphy  In  the  Sunwapta  Pass  area.  A  useful  map 
can  be  found  in  the  Guidebook  for  the  Fifth  Annual  Field 

•N 

Conference  (1955)  published  by  the  Alberta  Society  of 
Petroleum  Geologists.  Rocks  in  this  area  range  in  age 
from  Cambrian  to  Mississippian.  In  only  a  few  places  is 
the  folding  extreme  and  rarely  does  faulting  complicate 
the  picture. 


, 
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The  cliff  of  shaly  limestone  on  the  west  side  of 
Sunwapta  River  below  Mt.  Kitchener  has  been  identified 
as  Middle  Cambrian  by  Allan.  The  overlying  strata  of 
buff  calcareous  shales  may  be  Upper  Cambrian.  About 
ten  miles  to  the  southwest  near  the  Columbia  glacier 
P.  E.  Raymond  (1928)  described  two  Middle  Cambrian  trilo- 
bites.  The  strata  on  the  east  side  of  the  road  above 
the  Chalet  are  probably  Middle  Cambrian  shales. 

The  Athabaska  glacier  flows  across  the  axis  of  a 
long  trending  anticline.  With  minor  offsets  this 
anticline  runs  to  the  east  of  Mt,  Edith  Cavell,  45  miles 
to  the  northwest,  through  the  glacier,  and  continues 
south  of  Mt.  Sarbach  30  miles  to  the  southeast.  A 
thick  succession  of  Lower  Cambrian  quartzites  are  exposed 
at  Edith  Cavell,  while  Middle  and  Upper  Cambrian  is 
seen  at  Mt.  Kitchener;  Ordovician  at  Mt.  Erasmus  and 
Devonian  at  Mt.  Sarbach.  It  is  evident  that  the  axial 
trend  is  N  40°  W  to  N  50°  W,  parallel  to  that  of  the 
Rocky  Mountain  front  in  this  area  and  that  the  anticline 
plunges  to  the  southeast.  The  anticlinal  axis  crosses 
the  glacier  at  Line  B  just*11  below  the  first  step  and  the 
limbs  dip  away  from  the  axis  at  10°  to  15°. 

The  Castle  Mountain  Syncline  passes  through  Nigel 
Peak  exposing  Banff  shale  of  Mississippi an  age.  The 
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structure  under  the  Sunwapta  River  has  not  been  determined 
but  Allan  mentions  the  possibility  of  a  fault  at  the 
west  base  of  Wilcox  peak.  This  may  tie  in  with  the 
small  thrust  fault  placed  by  Hughes  where  Tangle  Creek 
nins  into  the  Sunwapta  River.  (See  Hughes  ’Section  GG’.) 

3.  THE  GRAVITY  MEASUREMENTS  AND  THEIR  REDUCTION 
3. 1  The  Field  Work 

The  observation  on  lines  A  to  F  were  made  with  the 
University  of  British  Columbia's  Worden  Meter  No.  35* 
According  to  the  manufacturer  this  meter  has  a  dial 
constant  of  0.41874  milligals  per  scale  division.  Later 
in  the  season  the  University  of  Alberta's  Worden  gravity 
meter  No.  XP  0  arrived  from  factory  servicing  and 
observations  were  made  with  it  on  Lines  B  to  H  inclusive. 
The  dial  constant  is  0.2607  milligals  per  scale  division 
as  determined  by  Dr.  G.  D.  Garland.  There  are,  therefore, 
two  independent  readings  on  85  of  the  127  stations 
occupied.  Wherever  possible  readings  with  meter  XP  0 
were  used  in  the  calculations  as  the  dial  constant  is 
more  accurately  determined,  the  meter  is  more  sensitive, 
and  it  has  a  lower  drift  rate.  Readings  were  taken  over 
a  sufficiently  wide  range  to  see  that  the  calibration  of 
the  two  meters  agreed  to  two  significant  figures. 
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A  base  station  was  set  up  at  C  1  behind  the  campsite. 
The  location  was  on  morainal  material  about  a  hundred 
feet  from  bedrock  and  measurements  indicated  that  it  was 
not  significantly  affected  by  ablation  or  horizontal 
movement  of  the  ice.  The  horizontal  velocity  amounted 
to  1  foot  per  month  at  this  station.  The  instrument  and 
tripod  must  be  accurately  leveled  before  a  reading  can 
be  taken.  Maintaining  a  level  base  on  the  glacier  is 
almost  impossible  unless  the  metal  tripod  legs  are 
insulated  from  the  ice.  Following  recommendations  of 
members  of  Operation  Hazen  to  Ellesmere  Island,  plasticene 
was  tried  and  found  to  provide  excellent  insulation. 

The  Worden  gravimeter  Is  an  astatic  instrument 

built  of  fused  quartz  and  housed  in  a  partial  vacuum. 

It  measures  small  variations  in  the  vertical  component 

of  gravitational  attraction  by  balancing  the  torque 

produced  by  a  mass  on  the  end  of  a  horizontal  arm  with 

the  restoring  force  of  a  spring.  The  astatic  element 

consists  of  a  spring  which  acts  in  the  same  direction 

as  the  force  of  gravity  and  nearly  cancels  the  elastic 

res coring  force.  The  unit  of  measurement  is  0.001 
2 

cm. /sec.  and  Is  termed  the  milligal. 

The  primary  base  station  was  on  the  Banff -Jasper 
highway.  The  absolute  value  of  the  acceleration  of  gravity 
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was  obtained  from  measurements  by  Garland  and  Tanner 
(1957)  at  their  station  No.  459  (Gatehouse).  From  this 
point  the  new  base  station  at  C  1  was  set  up  by  a  system 
of  looping.  That  is,  observations  were  made  at  station 
459>  next  at  a  new  station  which  we  can  call  A  and 
then  again  at  459*  Proceeding  to  a  new  loop  observations 
were  made  at  stations  A,  B  and  A.  All  closed  loops 
which  established  bases  were  completed  in  less  than  one 
hour  and  the  instrumental  drift  was  determined  and 
corrected  using  a  linear  drift  rate.  From  the  behavior 
of  the  drift  curves  between  sequential  loops  it  is 
estimated  that  the  base  at  C  1  is  correct  within  0.1 
milligal. 

In  carrying  out  the  survey  over  the  glacier  the 
instrument  was  back-packed  over  rough  glacier  ice  and 
was  subjected  to  shocks  which  made  the  drift  rate  more 
erratic.  Ties  were  made  at  base  stations  at  least  once 
every  two  hours.  From  a  study  of  the  drift  curves  and 
a  comparison  of  the  values  obtained  with  the  second 
meter  it  appears  that  all  readings  are  within  0.2 
railligals  except  for  6  stations  which  do  not  agree  by 
0.3  to  0.5  milligals.  This  agreement  is  quite  good 
considering  that  measurements  with  the  two  meters  were 
made  two  to  three  weeks  apart  and  it  is  not  certain  that 
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ablation  which  varied  between  0  and  3  feet  has  been 
exactly  corrected  for.  Since  depth  determinations  were 
made  on  individual  traverse  lines  Independently  of 
observations  made  elsewhere,  the  accuracy  of  each  line 
should  be  considered  separately.  On  this  basis  errors 
between  stations  all  on  one  line  do  not  exceed  0.1  milli- 
gals . 

3.2  Surveying 

All  the  stations  were  surveyed  by  plane  table  and 
alidade  and  the  ties  around  closed  loops  were  within 
one  foot.  Many  extra  points  were  surveyed  on  the 
mountain  slopes  adjacent  to  the  glacier  so  that  accurate 
terrain  corrections  could  be  made.  The  movement  stakes 
were  surveyed  independently  by  S.  Paterson  with  a 
Wild  T  2  (57767)  Theodolite  and  these  measurements  have 
recently  been  reduced.  The  results  are  in  agreement 
within  one  foot.  All  the  surveying  is  based  on  the 
Water  Resources  Branch  Bench  Mark  No.  1  which  has 
coordinate : 


Latitude:  52°  13'  12" 

Longitude:  117°  13'  30" 

Elevation:  6468.04  feet 

Additional  horizontal  control  was  obtained  from  Bench 
Mark  17  near  the  camp.  Its  latitude  is  52°  11’,  13"; 
the  longitude  is  117°  15'  08" . 
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Attempts  at  using  a  barometer  as  an  altimeter  were 
not  successful  as  the  temperature  and  pressure  varied 
from  the  sides  to  the  middle  of  the  glacier.  The  measure¬ 
ments  were  generally  accurate  within  10  feet  so  the 
barometer  was  used  to  fill  in  topographic  details 
within  250  feet  of  a  station.  Routine  readings  were 
made  in  camp  which  has  an  altitude  of  7340  feet.  The 
altimeter  reading  at  camp  varied  between  7078  and  7640 
with  the  reading  occasionally  rising  by  400  feet  within 
a  few  hours  when  a  storm  was  impending. 

To  estimate  the  amount  of  horizontal  movement  down 
the  glacier  two  markers  were  set  up  on  bedrock  on  either 
side  on  the  glacier.  The  Wild  T-2  Theodolite  was  aligned 
with  these  markers  and  the  holes  for  the  stakes  on  Line  C 
were  drilled  within  0.1  feet  of  this  line.  Forty-five 
days  later  the  instrument  was  reset  in  its  original 
position  and  the  original  line  marked  on  the  glacier 
opposite  each  stake.  The  down -dip  movement  of  the 
stakes  was  then  measured  simply  with  a  tape.  The  results 
are  shown  in  Table  I,  when  reduced  to  a  30  day  month. 

The  thickness  of  the  ice  along  the  flanks  was 
determined  from  this  gravity  survey  while  the  maximum 
depths  are'  from  a  borehole  200  feet  from  the  line.  There 
is  a  qualitative  relationship  between  the  depth  and  the 
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TABLE  I.  Surface  Movement  Studies  on  C  Line 


Stake 

Number 

Velocity 
(feet  per  month) 

Vertical 
of  Ice 

Thi ckness 
in  Feet 

C 

1 

1.0  +  0.1 

140  + 

20$ 

C 

2 

4.9  0.1 

330 

20$ 

C 

3 

10.0  0.1 

550 

20$ 

C 

4 

11.6  0.1 

700 

20$ 

c 

5 

12.6  0.1 

830  + 

10$ 

c 

6 

13-1  0.1 

940 

10$ 

c 

7 

13-7  0.1 

990 

10$ 

c 

8 

14.0  0.1 

1020 

10$ 

c 

9 

14.2  0.1 

1020 

10$ 

c 

10 

14.0  0.1 

980 

10$ 

c 

11 

14.3  0.1 

910 

10$ 

c 

12 

13.1  +  1.0 

740  + 

0 

OJ 

c 

13 

11.7  ±  0.1 

640 

20$ 

c 

14 

5.7  0.1 

420 

20$ 

c 

15 

0  0.1 

160 

20$ 
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surface  velocity  on  this  line.  A  major  problem  on  many 
glaciers  which  have  an  extensive  lateral  moraine  is  to 
determine  the  zero  Ice  thickness  line.  It  is  evident 
that  velocity  measurements  of  this  type  can  locate  this 
line  within  100  feet.  Difficulties  may  be  encountered 
in  the  lower  portions  of  the  glacier  where  stagnant  ice 
may  be  present.  Although  the  gravity  measurements  were 
made  on  a  moving  surface ,  this  study  indicates  that  the 
movement  was  not  enough  to  warrant  a  latitude  correction. 

The  rate  of  ablation  was  determined  from  measure¬ 
ments  on  the  movement  stakes  which  were  drilled  into  the 
ice  a  depth  of  7  feet.  The  ablation  during  July  varied 
from  0  to  5  feet  per  month,  the  amount  depending  on  the 
location  of  the  stake.  Since  the  surveying  and  the 
gravity  measurements  were  made  one  to  three  weeks  apart, 
a  correction  to  the  altitude  was  necessary. 

3*3  The  Reduction  of  the  Gravity  Measurements 

The  principal  facts  for  the  gravity  stations  are 
set  out  in  Table  II.  This  table  gives  the  station,  line 
number,  longitude,  latitude,  elevation,  observed  gravity 
and  the  various  types  of  anomaly.  It  is  necessary  to 
correct  for  variations  of  gravity  with  latitude,  eleva¬ 
tion  of  the  station  above  the  reference  surface  or  geold 
at  sea  level  and  the  attraction  of  the  surrounding 
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topography.  The  remaining  figure  is  termed  the  Bouguer 
anomaly  and  is  related  to  the  anomalous  mass  variations 
within  the  earth. 

The  variation  of  gravity  with  latitude  is  given  by 
the  International  Gravity  Formula  (1930). 

7  =  973.0490(1  +  0.0052384  sin2  0  -  0.0000059  sin2  2  0)  cm/sec2 
0  =  latitude  of  station 

The  variations  amount  to  0.024  milligals  per  second  of 
latitude  (101.6  feet)  at  the  glacier. 

The  free  air  correction  amounts  to  -  .094057  H  milli¬ 
gals  where  H  is  the  altitude  in  feet  above  sea  level. 

This  corrects  for  the  decrease  in  gravity  with  height 
above  sea  level.  The  Bouguer  correction  amounts  to 
0.034057  H  milligals  if  the  density  of  the  material 
between  the  surface  and  sea  level  is  2.67  grams  per  cubic 
centimetre.  This  corrects  for  the  attraction  of  an 
infinite  horizontal  slab  of  material  between  the  geoid 
and  the  earth's  surface.  A  terrain  correction,  T,  is 
then  made  to  account  for  the  topographic  variation  of 
the  surface  above  and  below  the  level  of  the  station. 

The  complete  expression  for  the  Bouguer  anomaly  is: 

g  -  ^Observed  “  Y  +  09*057  H  -  0.03*057  H  +  T 

The  expression  is  termed  a  simple  Bouguer  anomaly  if 
no  terrain  correction  is  made.  Figure  III  Is  a  contour 
map  of  the  Bouguer  anomalies.  (See  Appendix.) 
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To  correct  for  the  topography,  the  area  surrounding 
a  station  Is  divided  into  cylindrical  zones  bounded  by 
radii  r^  and  r2*  Each  annular  zone  is  subdivided  into 
N  sectors.  The  attraction  of  each  sector  is  given  by 
the  formula:  (Heiland,  p.  149). 


where  h  is  the  average  height  of  the  sector  above  or 


below  the  station,  S  is  the  density  and  k  is  the  gravi¬ 


tational  constant. 

Kayford's  values  for  the  radii  were  used  and  correc¬ 
tions  were  made  to  all  stations  for  his  zones  A  to  J 
inclusive.  This  corrected  for  the  topographic  effect  of 
the  mountains  and  valleys  to  a  distance  of  7*7  miles. 

The  effect  of  the  terrain  beyond  this  distance  is 
nearly  the  same  over  the  entire  glacier  and  does  not 
affect  the  accuracy  of  the  relative  gravity  anomalies. 
Bullard  (1935)  has  published  a  table  evaluating  the 
above  formula  for  different  values  of  r  and  h.  This 
table  is  not  accurate  enough  for  modern  gravity  surveys 
and  a  new  table  was  computed  for  use  here. 

Estimating  the  heights  in  each  of  the  9 6  sectors 
surrounding  a  station  involved  the  greatest  amount  of 


labour.  Care  had  to  be  taken  to  distinguish  between 
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sectors  containing  ice  which  has  a  density  of  0.90  and 
sectors  containing  rock  with  a  density  of  2.67-  The 
reduction  of  terrain  corrections  for  125  stations  required 
about  one  month  to  complete.  The  corrections  for  Kay ford 
zones  B  to  J  are  set  out  in  Table  III.  The  terrain 
corrections  were  very  large  because  of  the  proximity  of 
Kt.  Athahaska  (11,452' ),  Snow  Dome  (11,340'),  Mt.  Kitchener 
(11,500')  and  the  Sunwapta  valley  (6300').  The  maximum 
correction  was  14.9  milligals  and  the  minimum  was  6. 31 
milllgals . 

3.4  Density  of  the  Ice  and  the  Rock  Samples 

The  density  of  fresh  fallen  snow  is  about  0.06  to 
0.l6  gms.  per  cc.  and  by  a  complicated  process  of  melting 
and  refreezing  it  forms  neve'  with  a  density  of  0.45. 

Since  the  glacier  was  entirely  free  of  snow  or  neve 
when  the  gravity  measurements  were  made,  no  correction 
was  necessary  for  this  light  material. 

The  theoretical  density  of  pure  ice  can  be  calculated 
from  a  knowledge  of  the  lattice  constants  and  is  O.9168 
gm.  per  cc.  at  0°  C.  As  the  depth  of  burial  Is  increased 
the  air  bubbles  in  the  neve  are  expelled;  however,  some 
are  trapped  and  gaps  are  formed  between  the  ice  crystals 
so  that  the  theoretical  value  Is  never  realized. 

Mercantour  (1950)  has  studied  the  subject  and  concludes 
that  a  good  figure  to  use  for  glacier  Ice  Is  0.90  gm.  per  cc. 


Rock  samples  were  collected  from  the  area  on  either 
side  of  the  glacier  and  the  results  of  the  density 
measurements  are  given  in  Table  V.  The  samples  are 
listed  in  approximate  stratigraphic  order  with  the 
youngest  rocks  at  the  top  of  the  Table.  The  samples 
are  probably  all  of  Cambrian  age. 

The  mean  density  by  rock  type  is  as  follows: 

TABLE  IV.  Mean  Density  of  the  Rock  Types 

Number  of  Samples  Rock  Type  Mean  Density 

  gm/cc. 


15 

Limestone 

and 

Dolomite 

2.69 

9 

Shale 

2.79 

4 

Sandstone 

and 

Quartzite 

2.6l 

It  is  difficult  to  decide  if  this  suite  of  rocks  is 
representative  of  the  area  in  which  terrain  corrections 
were  made.  Limestone  is  the  predominant  rock  which 
outcrops  in  this  area  and  it  appears  that  the  density 
of  2.67  gm.  per  cc.  used  in  the  corrections  is  very 
nearly  the  correct  one. 

It  is  expected  that  Lower  Cambrian  quartzites 
appear  at  some  level  beneath  the  glacier  surface.  Incon¬ 
clusive  evidence  to  this  effect  was  obtained  from  a 
seismic  refraction  line  shot  along  the  glacier.  It  was 


' 


' 
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TABLE  V.  Densities  of  Rock  Samples 


Rock  Type 

Locality 

Density 

gm./cc. 

CAMBRIAN  SYSTEM 

Limestone,  black,  with  calcite  veins, 
cliff  forming 

Above  B  10 

2.69 

2.71 

Limestone,  black,  fine  grained 

Above  BM  17 

2.68 

Shale,  dark  grey,  calcareous 

Above  BM  17 

2.82 

2.83 

Limestone,  black,  fine  grained 

Basal  unit  on  south  side 

Above  3M  17 
Above  BM  17 

2.71 

2.72 

Limestone,  light  grey,  massive,  hard 

Opposite  3  10 

2.69 

Dolomite,  dark  green,  hard,  well  bedded 
about  300'  thick 

Above  M  23 

2.63 

2.67 

Limestone,  light  grey,  and  chexty 

dolomite, 
weathers  yellow 

Opposite  M  14 

2.82 

2.83 

Limestone,  b^ack,  massive,  with  calcite 
veins.  About  50  feet  thick 

Above  D  L  on 
north  side 

2.73 

2.69 

Limestone,  dark  grey,  shaly 

Above  D  L  on 
north  side 

2.60 

2.6l 

Shale,  grey,  well  bedded 

Above  D  L  on 
north  side 

2.72 

Limestone,  black,  poorly  preserved 
trilobite  fragments. 

Above  D  L  on 
north  side 

2.63 

Sandstone,  quartzite,  well  cemented, 
light  yellowish  grey 

Above  D  L  on 
north  side 

2.58 

2.63 

Shale,  dark  grey 

Above  D  L  on 
north  side 

2.76 

2.79 

2.81 

2.77 

Shale,  dark  grey 

Above  D  L  on 
north  side 

2.83 

Quartzite,  pinkish  white 

Basal  unit  on  south  side 

Above  D  L  on 
north  side 

2 . 66 
2.58 

MEAN  DENSITY 

2.71 

, 

• 
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hoped  that  arrivals  from  a  high  velocity  layer  such  as 
a  bed  of  Middle  Cambrian  limestone  would  be  detected  by 
the  refraction  program.  'Tentative  study  of  the  time- 
dlstance  plot  indicates  that  first  arrivals  travelled 
with  a  velocity  of  12,200'/s  through  the  ice  while  the 
secondary  arrivals  indicate  velocities  no  higher  than 
14,000  feet  per  second.  This  is  a  velocity  that  would 
be  expected  from  a  layer  of  sandstone  but  is  far  too 
loiv  for  a  bed  of  consolidated  limestone.  The  mean 
density  of  the  sandstones  is  lower  than  was  expected 
and  there  may  have  been  some  leaching  which  Increased 
the  porosity.  However,  the  density  is  only  2%  lower 
than  the  value  of  2.67  which  was  used  in  calculating 
the  thickness  of  the  ice  column.  The  density  difference 
(between  rock  and  ice)  used  in  all  computation  was 
1.77  gm.  per  cc*  and  this  figure  appears  to  be  correct 
within  2%. 

3*  5  Accuracy  of  the  Data 

The  relative  field  measurements  on  each  line  are 
accurate  to  within  0.1  milligals.  The  elevations  are 
known  to  a  precision  of  +  1  foot  so  that  the  combined 
free  air  and  Bouguer  corrections  are  accurate  to  +  0.06 
milligals.  The  principal  error  and  the  one  that  is 
most  difficult  to  evaluate  is  that  due  to  the  terrain 
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corrections.  The  uncertainty  in  the  density  contributes 
an  error  of  about  2$.  This  produces  an  error  of  0.1 
milligals  between  a  station  in  the  middle  and  a  station 
at  the  sides  of  the  glacier.  For  each  station  there 
are  9 6  sectors  each  of  which  is  subject  to  an  error  of 
about  0.01  milligals.  These  reading  errors  are  in  the 
nature  of  a  "random  walk".  In  the  measured  quantity  T 
there  are  96  causes  of  error,,  each  of  which  will  add  the 
amount  0.01  or  -0.01  with  equal  probability.  The 
probability  of  an  error  of  +  0.04  is  only  0.2  while  the 
probability  of  an  error  as  large  as  +  0.1  milligals  is 
0.004.  • 

In  the  preceding  discussion  it  has  been  tacitly 
assumed  that  the  topographic  maps  are  perfect.  In  some 
of  the  zones  a  determinate  error  may  be  introduced  due 
to  a  systematic  distortion  of  the  contours.  A  large 
scale  map  on  the  scale  of  1:4800  was  constructed  and 
contoured  in  the  field.  Zones  B  and  C  which  cover  an 
area  with  a  radius  of  755  feet  were  entirely  corrected 
using  this  map  and  it  is  unlikely  that  any  significant 
error  was  introduced  here.  Zones  F  to  J  which  cover 
an  annular  area  1.4  to  7*7  miles  from  a  station  were 
also  corrected  on  maps  of  adequate  scale  and  accuracy. 
The  1919  topographic  maps  on  a  scale  of  1:62,500  are  not 
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adequate  for  correcting  zones  D,  E  and  P.  The  magnitude 
of  the  error  can  only  be  roughly  determined  from  plausible 
changes  in  the  shape  of  the  contours.  It  is  estimated 
that  the  error  is  no  more  than  +0.2  milligals.  The 
terrain  corrections  are,  therefore,  thought  to  have  a 
precision  of  at  least  +  0.4  milligals  and  the  precision 
of  the  Bouguer  anomaly  is  +0.5  milligals. 

4.  THE  INTERPRETATION  OF  GRAVITY  ANOMALIES  OVER  A  GLACIER 

4. 1  Introduction 

The  interpretation  of  gravity  anomalies  involves 
an  evaluation  of  the  mass  distribution  causing  the 
anomaly.  A  unique  interpretation  cannot  be  found  as 
this  is  a  problem  in  potential  theory,  and  one  can 
rearrange  the  subsurface  mass  distribution  in  an  infinite 
number  of  ways  to  reproduce  the  same  gravitational  field. 
In  the  present  problem  there  is  additional  information 
which  places  restriction  on  the  interpretations  possible. 
The  density  difference  producing  the  anomaly  is  known 
within  2%.  If  the  shape  of  the  bedrock  profile  were 
known,  the  central  depth  to  bedrock  would  be  the  only 
unknown  variable  and  it  could  be  determined  exactly 
from  a  knowledge  of  the  gravitational  field.  Sometimes 
the  depth  of  a  glacier  will  be  known  from  a  borehole  or 
from  seismic  work.  In  this  case  the  shape  of  the  bedrock 
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profile  can  be  determined  approximately.  The  detection 
of  small  sinusoidal  variations  on  the  bedrock  would  be 
limited  by  the  accuracy  of  the  gravity  data. 

Determining  the  thickness  of  a  glacier  with  gravity 
data  was  first  attempted  in  1955  on  the  Austerdalsbre, 
a  Norwegian  glacier,  by  C.  Bull  and  J.  R.  Hardy.  No 
topographic  maps  were  available  so  conventional  terrain 
corrections  could  not  be  made.  An  estimate  of  the  effect 
was  made  by  three  traverses  at  the  foot  of  the  glacier. 
The  thickness  was  calculated  by  assuming  that  the  sheet 
of  ice  extended  to  infinity  in  a  horizontal  direction. 
Essentially  this  is  a  Bouguer  correction  which  is  useful 
in  mapping  an  ice -sheet  as  shown  by  C.  A.  Littlewood 
on  Baffin  Island  In  1952.  The  method  can  seldom  be 
applied  to  a  glacier  as  it  is  a  three  dimensional,  or 
at  best,  a  two  dimensional  body.  The  Austerdalsbre 
could  be  regarded  as  an  ice  sheet  of  infinite  extent 
because  it  was  less  than  100  meters  thick  and  about 
1000  meters  wide.  The  principal  errors  arise  from  the 
method  of  making  terrain  corrections.  The  authors 
claim  that  these  are  correct  within  1  to  2  milligals  in 
the  center.  Errors  in  the  ice  thickness  are  given  as 
20$  on  the  lower  two  profiles  and  40$  on  the  upper  two. 

Thiel,  LaChapelle  and  Behrendt  made  gravity  measure¬ 
ments  on  the  Lemon  Creek  Glacier,  Alaska  In  1956.  No 
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topographic  maps  were  available  and  terrain  corrections 
were  made  approximately  with  a  two  dimensional  model  of 
the  mountain,  the  model  extending  to  infinity  in  the 
third  direction.  A  correction  was  made  for  the  regional 
trend  by  comparing  the  Bouguer  anomaly  of  the  end 
stations.  As  at  the  Austerdalsbre,  a  first  approximation 
was  made  by  applying  the  infinite  slab  Bouguer  correction. 
However,  in  this  case,  a  second  approximation  was  made 
by  a  line  integral  method  of  calculating  the  gravity  of 
any  desired  cross  section  (Hubert,  1948).  The  Lemon 
Creek  glacier  proved  to  be  a  little  over  200  meters 
deep  but  the  authors  do  not  give  an  estimate  of  the 
accuracy. 

Jacobs,  Grant  and  Russell  made  gravity  measurements 
on  the  Salmon  Glacier,  3.C.  This  is  a  deep  glacier 
with  a  borehole  reaching  a  depth  of  725  meters.  The 
interpretative  technique  was  a  two  dimensional  variant 
of  the  one  used  by  Thiel,  et  al.  The  accuracy  of  the 
data  is  discussed  fully  and  the  thickness  is  believed  to 
be  correct  within  -10$  and  4-25$.  Most  important  was  a 
comparison  with  seismic  and  borehole  data  on  one  line. 

The  agreement  between  these  data  was  within  10$. 

4.2  Theory  of  Gravity  Calculations 

The  potential  at  the  origin,  P,  due  to  a  mass  of 
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density,  §  ,  distributed  evenly  throughout  a  volume,  T  , 

not  occupying  the  origin  is : 


u 


kS 


dx  dy  dz 
'x?+  y2  +  z 


2 


(1) 


A  left  handed  cartesian  coordinate  system  is  employed, 
with  the  z  axis  vertically  downward.  The  gravitational 
constant  is  denoted  by  the  symbol  "k".  If  we  let  r^  =  x^ 
and  change  to  cylindrical  coordinates  (r,  ip,  z),  (l) 
becomes : 
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The  gravity  anomaly  due  to  a  pie -shaped  prism  of 
thickness  dz  and  radius  r  is° 


Figure  14 


The  gravity  anomaly  due  to  a  lamina  of  infinitesmal 
thickness  dz  is: 
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Ag  =  V  dz 


(3) 


where 


V  =  kS  {  |  dij/ 


di//}  W 


If  the  radius  of  the  lamina  becomes  very  large  the 


second  integral  approaches  zero.  Then, 


Ag  -  2  7AkS  dz 


(5) 


which  is  the  Bouguer  correction  formula  for  an  Infinite 
bed  of  thickness  dz.  Formulas  (3)  and  (4)  were  programmed 
for  the  Royal  McBee  LGP-30  digital  computer  to  obtain 
the  gravitational  attraction  of  the  glacier  at  any 
station  following  the  method  developed  by  Talwani  and 
Ewing  (i960).  These  formulas  are  also  the  analytic 
basis  of  a  mechanical  integrator  for  the  computation  of 
gravity  anomalies  (Siegert,  1942). 

Consider  next  the  attraction  due  to  a  line  element 
extending  to  plus  and  minus  Infinity  in  the  y  direction. 
From  ( 1 ) 
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(6) 


y  =A/  x2  +  z2  tan  Q; 
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Let  r2  =  x2  4-  z2 


dy  -  /  X2  4-  ,2  sec2  O  d© 


V  x  +  z 


Ag 


2kSz  f 
.2^2  / 


x2+  z2  J 


cos  6?  d@ 


(7) 
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From  this  formula  it  is  easy  to  calculate  the  effect 
of  a  two  dimensional  body  such  as  a  rectangle  (Heiland, 
p.  150).  We  can  then  build  up  a  two  dimensional  model 
of  a  glacier.  The  gravitational  attraction  of  the  shaded 
rectangle  is  given  by: 


b  D 
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Figure  15 
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If  (7)  is  converted  to  polar  coordinates  with 


Ag  =  2kS  —  r  dr  d<t)  (9) 


=  2k8(rn„-'-n)(cos  t,r  cos  4>J 


In  graphical  interpretation  this  formula  is  used  to 
construct  a  graticule  chart.  The  gravity  anomaly  is 
obtained  by  superimposing  the  chart  on  a  cross'  sectlon 
of  the  glacier  with  the  origin  at  the  gravity  station 
and  counting  the  sectors  within  the  outline  of  the 
glacier.  Nettleton  (p.  117)  provides  a  useful  formula 
and  table  to  give  an  "end  correction"  for  reducing  the 
gravity  effect  of  a  two  dimensional  body  to  a  body  of 
finite  length. 

4.3  A  Two  Dimensional  Model  of  the  Glacier 

The  process  of  finding  the  depth  and  shape  of  a 
glacier  is  essentially  a  trial  and  error  process. 

Formula  (8)  for  a  two  dimensional  model  with  rectangular 
cross  section  can  be  adapted  as  shown  in  the  diagram 
below  the  formula  to  yield  the  gravity  anomaljr  of  any 
arbitrary  shape  to  any  desired  accuracy  by  choosing 
(D  -  d)  small  enough.  Alternately  a  graticule  chart 
can  be  made  using  formula  (9)  with  (r  ^  -  rn)  *  1  unit 
and  (cos  -  cos  0n)  =  0.05*  The  gravity  anomaly  of 
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any  two  dimensional  body  of  arbitrary  cross  section  can 
be  found  in  3  to  5  minutes  per  station.  This  is  con¬ 
siderably  faster  than  computation  with  formula  (8). 

Geomorphologists  have  described  glacier  valleys  as 
being  U-shaped  or  parabolic.  H.  Svensson  (1959)  has 
demonstrated  that  the  glacial  valley  Lapporten  in 
northern  Sweden  can  be  described  by  an  equation  of  the 
form 


=  Axh 


y 


in  which  b  is  very  nearly  2.  It  was*  therefore,  decided 
that  a  good  first  approximation  to  the  shape  of  the 
bedrock  beneath  the  Athabaska  glacier  would  be  a 
parabola  of  the  form: 


or  an  ellipse: 


The  origin  is  on  the  surface  in  the  centre  of  the  glacier. 
The  semi -major  axis  is  given  in  terms  of  "a",  the 
distance  from  the  centre  to  the  edge  of  the  glacier,  and 
the  semi -minor  axis  is  ’'bn,  the  central  depth.  These 
two  models  proved  highly  useful,  particularly  as  the 
parabola  gives  a  small  side -wall  angle  while  the  ellipse 
provided  a  model  with  very  steep  sidewalls. 

Figure  5  gives  a  comparison  of  the  shape  and  theor¬ 
etical  gravity  anomaly  of  a  rectangular,  elliptical  and 
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parabolic  cross  section  for  the  glacier  at  the  approximate 
position  of  the  borehole.  One  of  the  points  to  notice 
is  that  the  gravity  anomaly  is  still  4  to  6  milligals 
where  the  glacier  has  zero  thickness.  The  zero-thick¬ 
ness  anomaly  varies  considerably  with  both  depth  and 
shape.  The  zero  thickness  line  can  seldom  be  determined 
very  accurately  and  it  is  often  so  close  to  the  steep 
mountain  sides  that  a  gravity  station  cannot  be  located 
on  it  because  the  terrain  correction  will  be  uncertain 
by  several  milligals.  Figure  5  also  indicates  that  the 
gravity  anomaly  has  the  greatest  rate  of  change  near 
the  edge  of  the  glacier.  These  considerations  will 
usually  make  it  unwise  to  tie  the  Bouguer  anomaly  from 
the  field  data  to  the  end  points.  The  Bouguer  anomaly 
does  not  give  us  information  about  the  value  of  the 
zero  ordinate  for  the  glacier  anomaly.  All  we  have  Is 
a  relative  anomaly  from  the  centre  to  the  end  station. 

The  centre  points  have  been  tied  together  in  any 
comparison  of  gravity  anomalies  from  field  or  model 
studies.  Figure  5  indicates  that  if  the  depth  to  the 
centre  is  known,  the  gravity  data  can  predict  the  shape 
with  considerable  success. 

The  choice  of  depth  on  each  line  was  guided  by  the 
information  from  one  borehole  and  from  a  knowledge  of 
the  seismic  results  (personal  communication  from  Professor 
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Figure  5*  A  comparison  of  the  theoretical  gr  vity  anomalies  of 
•»  three  two-dimensional  bodies*  The  bodies  are  4200  feet  wide  along 

the  surface  (x  direction)}  1024  feet  deep  along  the  z  direction  and 
extend  to  infinity  in  the  y  dlredtion*  The  density  contra  t  is  1*77  gu/cc* 
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J.  S.  Savage).  However,  various  other  depths  were 
modelled  to  check  the  seismic  results  and  to  obtain  a 
value  for  the  precision  with  which  depth  can  be  predicted 
by  gravity  results  alone.  Based  on  this  work/  a  shallower 
depth  was  chosen  on  Lines  A,  E.  H  and  G.  The  shape  of 
the  cross  section  was  then  determined  on  the  basis  of 
the  gravity  depth  rather  than  the  depth  indicated  by  a 
preliminary  Interpretation  of  the  seismic  work. 

As  shown  on  page  32,  the  theoretical  gravity  anomaly 
of  a  two  dimensional  body  with  a  parabolic  cross  section 
can  be  approximated  as  closely  as  required  by  the  repeated 
use  of  the  formula  for  the  gravity  of  a  rectangular 
prism  (formula  8).  This  expression  can  easily  be 
programmed  for  a  digital  computer  to  yield  the  theoretical 
gravity  at  any  given  point  along  a  line  on  the  surface. 

Let  the  z  axis  pass  through  the  center  of  the  parabola 
in  Figure  15.  The  symbols  used  will  conform  with  those 
in  formula  8  except  for  the  appropriate  subscripts. 

h  =  maximum  depth  of  the  glacier 

a  =  half  width  of  the  glacier  at  the  surface 

t  =  Dj_  -  d^,  the  thickness  of  a  rectangular  prism 

=  coordinate  of  observation  point  from  the  center 
of  the  glacier 

n  =  number  of  observation  points 
i  =  h/t,  number  of  rectangular  prisms 
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The  equation  of  the  parabola  is 

z  =  Ax2  -f  C 
h  2 

or  z  =  -  — 5—  x^  -i-  h 

a^ 

L  =  x  =  a,^  1  _  zi 
2  ^  IT 

D±  =  Zf  +  t 

di  =  Zi  -  t 

For  a  particular  observation  point,  £>n  \ 

xi  =  I  £„ I +  | 

The  program  is  designed  to  choose  a  given  depth,  h,  and 
observation  point,  Cn  >  and  compute  formula  (8)  for  each 
of  the  i  rectangular  prisms.  It  sums  these  up  and  prints 
out  the  answer.  The  process  is  repeated  for  all  n 
observation  points  and  for  each  maximum  depth,  h. 

A  program  involving  the  above  equation  was  written 
for  the  LGP-30  digital  computer  and  carried  out  for  15 
observation  points  along  Line  C.  The  thickness  of  the 
prisms  was  chosen  to  be  50  feet.  For  a  determination  of 
the  best  fit,  the  Bouguer  anomaly  of  the  center  station 
was  tied  to  the  theoretical  gravity  value.  The  root 
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mean  square  of  the  deviations  from  the  theoretical  values 
was  computed  for  each  of  the  maximum  depths.  The  results 
are  summarized  in  Table  VIII. 

TABLE  VIII  Results  of  Two  Dimensional  Computer  Program 

Line  C  -  15  observation  points 

Depth  Root  Mean  Square  of  Deviations  of 

Theoretical  from  Observed  Gravity 
(feet)  (milligals) 


800 

0.80 

900 

0.67 

1000 

0. 54 

1100 

0.  51 

1200 

0.  54 

1300 

o.6o 

From  a  seismic  program  and  from  two  nearby  boreholes 
which  reached  bedrock,  the  depth  is  known  to  be  1050 
+  30  feet.  On  the  basis  of  the  above  analysis  a  depth 
of  1100  would  have  been  chosen. 

The  procedure  illustrated  above  can  be  used  in  a 
systematic  way  on  the  central  part  of  any  glacier  for 
which  the  observed  gravity  anomaly  is  known  to  the 
required  accuracy.  The  expression  for  the  gravity  of  a 
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rectangular  prism  can  be  adapted  to  an  elliptical  cross 
section  with  equal  ease.  If  the  Bouguer  anomalies  are 
asymmetrical,  a  similar  degree  of  asymmetry  can  be 
programmed  into  the  basic  parabolic  or  elliptic  cross 
section  and  the  maximum  depth  determined  from  the 
smallest  root  mean  square  of  the  deviations  as  in  the 
example  given  above. 

4.4  Computer  Program  for  the  Gravitational  Attraction 

of  Body  With  Finite  Dimensions 

Once  a  two  dimensional  approximation  has  been  found 
for  the  glacier  the  body  of  ice  can  be  treated  as  an 
irregular, three  dimensional  ore  body  whose  gravity 
anomaly  can  be  computed  by  a  method  developed  by  Taiwan! 
and  Ewing.  (See  the  top  diagram  in  Figure  6.)  A 
contour  map  of  the  model  of  the  glacier  is  drawn  up  and 
each  contour  is  replaced  by  a  horizontal  irregular  n-sided 
polygonal  lamina  as  shown  in  the  middle  diagram  of 
Figure  6.  The,  contour  interval  used  was  100  feet.  The 
polygon  should  fit  the  contour  map  closely  if  an  observa¬ 
tion  point  is  nearby.  Let  m  be  the  number  of  observa¬ 
tion  points  with  coordinates  (  £,  7?j  )  and  elevation  Ej, 

1  j  ^  m.  Let  n  be  the  number  of  points  which 

define  the  comers  of  the  lamina  with  coordinates  (X^,  Yj_ ) 
and  elevation  Ej/  1  -  i  n.  The  coordinates  of  the 
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Figure  6.  Geometry  for  the  evaluation  of  the  gravity  anomaly 
caused  by  the  density  contrast  between  ice  and  bodbrock.  V  is 
gr  vity  ano  laly  at  a  point  (A-5)  due  to  a  horizontal  lamina  at 


depth  z 
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lamina  are  translated  so  the  origin  of  the  cartesian 
axis  is  at  the  observation  point  (£.  rj.  ). 

xi  «  Xi 

y±  ~  Y i  - 

zi  *  Ej  -  Ei 

Following  Taiwan!  and  Ewing,  formula  (4)  was  reduced 
to  a  form  suitable  for  evaluation  by  a  digital  computer. 
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1  if  is  positive. 


3 


-  1  if 


is  negative. 


W  «  ♦  1  if  is  positive.  W  «  -  1  if  is  negative. 


The  program  prints  out  a  table  of  ZJ  ^  iV,  i  a™1 
V  for  each  lamina  and  then  reads  the  data  for  the 
next  lamina. 

A  glacier  is  a  unique  case  of  an  "ore1  body  because 
the  observation  points  straddle  the  "ore”  body  itself. 
Therefore,  only  the  lamina  below  the  observation  point 
need  be  computed.  The  material  above  the  level  of  the 
observation  point  has  been  removed  by  the  terrain 
corrections. 

As  written  for  the  LGP-30  computer,  the  program  is 
in  fixed  point  and  the  coordinates  of  the  observation 
points  should  be  less  than  30,000  feet.  The  time  to 
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compute  one  lamina  is  10 n  seconds  per  observation  point, 
where  n  is  the  number  of  sides  to  the  lamina.  The  value 
of  kS  is  6.667  •  10"^  cm-Vgm.  sec^.  1.77  gm/crn^  or 
3.597  •  10-3  milligals  per  foot. 

Two  models  were  computed  for  each  of  the  8  lines. 

For  the  majority  of  the  lines,  the  second  model  was  20% 
shallower  up-dip  from  the  line  in  question.  In  addition, 
2  extra  models  were  determined  for  lines  3  and  D  to 
provide  a  better  fit.  The  total  computer  time  for  45 
observation  points  for  the  two  models  was  13-1/2  hours. 
Two  hours  of  extra  computer  time  was  required  to  adjust 
the  fit  on  lines  B  and  D.  Including  the  time  required 
to  determine  the  coordinates  of  the  observation  points 
and  lamina,  the  program  was  not  completed  significantly 
faster  than  it  could  have  been  done  with  Siegert's 
mechanical  integrator.  However,  modification  to  the 
model  requires  little  change  to  the  input  data  and  can 
be  done  very  quickly.  The  LGP-30  is  a  very  slow  computer 
and  a  larger  model  would  complete  the  task  in  a  fraction 
of  the  time. 

Formula  (3)  for  the  gravity  anomaly  caused  by  the 
whole  glacier  was  not  programmed.  The  bottom  diagram  of 
Figure  6  illustrates  the  numerical  integration  for  one 
observation  point. 
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4. 5  Results 

The  results  of  the  computer  program  are  presented 
in  Table  VI  and  Figures  7  to  12  inclusive.  The  top 
diagram  is  a  transverse  cross  section  of  the  glacier 
together  with  various  other  models  computed  for  the 
purposes  of  comparison.  The  lower  two  graphs  show  a 
comparison  of  the  observed  gravity  anomaly  and  the 
theoretical  gravity  anomaly  of  the  model.  The  observed 
Bouguer  gravity  values  did  not  have  to  be  corrected  for 
a  regional  gradient  because  the  cross  sections  are  along 
the  strike  of  geologic  structures.  This  is  confimed  by 
a  regional  gravity  survey  in  the  area  (Garland  and 
Tanner,  1957).  The  final  model,  shaded  in  grey,  was 
generally  interpolated  from  the  computed  model  to  give 
a  better  fit  than  is  illustrated  by  the  anomaly  curve. 

The  A  line  appears  to  be  slightly  shallower  than 
the  seismic  data  would  indicate.  However,  the  seismic 
shot  point  was  off  the  line  at  L  10  so  it  is  not  possible 
to  make  an  exact  quantitative  comparison.  Transverse 
cross  sections  were  determined  seismically  on  Lines  D  and 
C  and  this  data  agrees  within  +  100  feet.  The  main 
difference  is  that  the  cross  section  determined  by 
gravity  is  notably  smoother  than  the  seismic  section. 
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Figure  7.  Cross  section  of  the  Athabaska  Glacier  together  with  the 
theoretical  and  observed  gravity  anomaly.  Model  Ill^ia  the  best  fit 
according  to  gravity  data.  Horizontal  and  Vertical  scales  are  equal 
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Figure  8.  Cross  section  of  the  Athabaslca  Glacier  together  with  the 
theoretical  and  observed  gravity  anomaly.  Model  III  M  is  the  best 
fit  according  to  gravity  data.  Horizontal  and  Vertical  scales  are 
equal. 
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Figure  9*  Gross  section  of  the  Athabaslca  Glacier  together  with  the 
theoretical  and  observed  gravity  anomaly.  Model  I  is  the  best  fit 
according  to  gravity  data.  Horizontal  and  vertical  scales  are  equal 
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Figure  10.  Cross  section  of  the  Athabasko  Glacier  together  with  the 
theoretical  and  observed  gravity  anomaly.  Model  III  is  the  best  fit 
According  to  gravity  data.  Horizontal  and  vertical  scales  are  equal. 
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Lines  E  and  H  are  10$  shallower  than  seismic  data 
indicates.  The  longitudinal  seismic  section  Indicates 
that  the  bedrock  is  quite  rough  in  this  section  and,  as 
expected,  gravity  data  tends  to  smooth  out  the  irregulari¬ 
ties.  The  longitudinal  profile  based  on  gravity  data 
(Figure  13)  is  everywhere  within  +  50  feet  of  the 
seismic  data  between  lines  B  to  D  and  within  +  100  feet 
between  lines  D  to  P.  Line  F  is  remarkable  for  its 
shallowness.  The  deeper  model  (450’)  which  was  programmed 
to  check  this  does  not  fit  the  field  data  but  the  fit  is 
quite  good  for  the  model  with  a  depth  of  370  feet 
(Figure  12). 

The  value  of  three  dimensional  computations  is 
best  illustrated  on  lines  near  the  terminus  of  the 
glacier.  It  was  possible  to  assign  an  absolute  value 
to  the  anomaly  produced  by  the  field  data  by  means  of 
a  station  at  the  foot  of  the  glacier.  The  terrain 
correction  here  is  certainly  accurate  within  0.2  milli- 
gals  and  the  absolute  depth  should  be  accurate  within 
+  10$.  The  seismic  record  has  a  poor  reflection  at  a 
depth  of  about  490  feet.  As  illustrated  on  the  bottom 
of  Figure  12,  this  is  incompatible  with  the  gravity 
results.  Models  were  programmed  for  a  depth  of  1/2  and 
1/3  the  seismic  depth  to  check  on  multiple  reflections. 

On  the  basis  of  the  gravity  study  the  reflection  appearing 
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Figure  11.  Cross  section  of  the  Athabaska  Glacier  together  with  the 
theoretical  and  observed  gravity  anomaly.  Model  III  is  the  best  fit 
according  to  gravity  data.  Horizontal  and  vertical  scales  are  equal. 


Figure  12.  Cross  section  of  the  Athabaska  Glacier  together  with  the 
theoretical  and  observed  gravity  anomaly.  Model  II  fits  best  for  the 
F  Line  and  Model  IV  for  the  G  Line.  Horizontal  and  Vertical  scales 
are  equal. 
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on  the  record  is  the  second  multiple.  The  primary 
reflection  should  appear  at  0.27  seconds  (+phase  lag  due 
to  filtering).  The  primary  reflection  and  its  first 
multiple  are  probably  obscurred  by  surface  waves. 

The  deepest  part  of  the  glacier  is  on  a  line  through 
B  6,  C  8  and  D  6.  This  is  best  seen  by  a  comparison  of 
the  Bouguer  Map,  Figure  3(Appendix),  and  the  Isopach 
Map,  Figure  4  (Appendix).  The  position  of  Line  L  was 
chosen  on  the  basis  of  a  week’s  study  of  the  stake 
movement  on  Line  C  (See  Table  l)  and  is  only  about  200 
feet  from  the  deepest  part  of  the  glacier. 

The  cress  section  of  the  glacier  at  lines  C  and 
E  is  very  nearly  a  parabola.  Line  F  has  a  cress  section 
midway  between  an  ellipse  and  a  parabola.  At  Line  H 
the  glacier  narrows  appreciably  and  the  shape  of  the 
cross  section  closely  resembles  an  ellipse.  Lines  B  and 
D  have  cross  sections  which  are  more  irregular. 

The  top  diagram  in  Figure  13  (Appendix)  shows  the 
Bouguer  anomaly  along  Line  L,  the  central,  longitudinal 
traverse.  It  gives  a  qualitative  representation  of  the 
thickness  of  the  glacier.  It  Is  now  possible  to  subtract 
the  gravity  anomaly  of  the  glacier  from  the  Bouguer 
anomaly  and  obtain  a  residual  Bouguer  anomaly  which 
should  reflect  the  underlying  geology.  (Use  was  made  of 
Table  VI  -  Results  of  Computer  Program,  Appendix) .  If 
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the  glacier  computations  have  been  properly  made,  the 
remaining  curve  should  be  very  smooth  and  become  more 
negative  to  the  southwest.  This  negative  trend  Into  the 
Rocky  Mountain  Trench  is  shown  on  a  Rouguer  Gravity  Map 
of  British  Columbia  and  Alberta  (Garland  and  Tanner, 

1957).  The  residual  anomaly  is  indeed  as  described  above. 
However,  between  lines  H  and  G  the  anomaly  decreases 
faster  than  the  regional  rate  and  indicates  less  dense 
material  under  the  bulk  of  the  glacier.  The  inflection 
point  is  on  Line  F  which  has  a  pronounced  bump  on  the 
glacier  floor.  These  features  may  be  related  to  the 
structural  geology  of  the  underlying  section  but  not 
enough  is  known  about  the  density  and  the  stratigraphic 
sequence  to  make  a  more  definite  statement. 

5.  APPLICATION  OF  THE  RESULTS  TO  GLACIOLOGY 

In  19^-9  Orowan  made  some  calculations  on  the  yield 
stress  of  the  ice  by  assuming  the  ice  behaved  as  an 
Idealized  plastic  solid.  J.  F.  Nye  (1951,  1957)  has 
developed  this  theory  much  further  and  derived  equations 
for  the  stress  distribution  and  the  velocity  of  a  two 
dimensional  slab  of  ice  flowing  down  a  gently  undulating 
rough  slope.  Deformation  is  assumed  to  be  negligible 
below  a  critical  yield  stress  and  the  shear  stress  does 
not  rise  above  this  value. 


■ 
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Allow  a  two  dimensional  wedge  of  ice,  with  surface 
slope,  a  ,  to  rest  on  a  rough  inclined  surface  of  slope,  p, 
as  shown  in  Ihe  diagram  below.  The  ice  cannot  slide  down 
except  by  internal  plastic  shear. 


Figure  16 


Following  Nye's  (1952  C)  analysis,  the  approximate  mean 
normal  pressure  on  AB  is  yOg  h  where  p  is  the 
density  and  g  is  the  gravitational  acceleration.  The 
normal  force  on  AB  is  \p  gh2  .  The  normal  force  on  DC 
is  X-p  gh2+.  -^(yOg^ySx.  The  force  parallel  to  the  bed 
is  77  (Sx)  where  r  is  the  shear  force  per  unit  area  exerted 
by  the  bed  on  the  ice  and  the  component  of  weight  of 
the  section  is  h(Sx)yOg  sin  /3  .  To  a  first  approxima¬ 

tion: 

^(-^/9gh2)8x  +  pq  h-/3'8x  =  T-Sx 


r 

p  gh 


or 


dh 
d  x 


+  (3 
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(^-  +  (3  )  is  the  slope  of  the  upper  surface,  CL  , 

therefore : 


T 


/OghQ 


(11) 


The  thickness  must  satisfy  the  following  simple  equation: 


where 


(12) 


Nye  (1952  C)  has  modified  this  formula  to  take  into 
account  the  sides  of  the  glacier  valley.  In  a  glacier 
of  arbitrary  but  constant  cross  section,  the  average 
value  of  shear  stress,  Ta  ,  at  the  bed  was  found  to  be: 


Ta 

R 


p  g  R  sin  CL 

A_ 

P 


(13) 


where  A  is  the  area  of  the  cross  section  perpendicular 
to  the  bed,  and  P  is  the  perimeter  of  this  cross  section. 
The  formula  is  approximately  true  if  the  cross  section 
varies  slowly  and  the  current  value  of  R  and  the  surface 
slope,  CL  ,are  used. 

If  the  cross  section  is  a  parabola: 
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If  the  cross  section  is  an  ellipse: 

P  =  TT 

A  =  ^  bh 

If  the  cross  section  is  irregular  in  shape,  approxi¬ 
mate  methods  must  be  used  to  find  the  perimeter  and 
area.  All  of  the  quantities  above  are  known  and  it  was 
possible  to  calculate  the  shear  stress  at  the  bed  on 
Lines  B,  C,  D,  E.  H  and  P.  The  gravitational  constant 

p 

was  taken  to  be  981  cm/sec  and  the  density  was  0.90 
gm/cm^ . 


TABLE  VII .  Average  Shear  Stress  on  the  Bedrock 


Line 

Depth  (h) 
(meters) 

Width  (d) 
(meters) 

R 

(meters) 

a 

0  1 

A  Ta 

(10  dynes/cm^) 

B 

326 

1055 

170 

3 

24 

0.9 

C 

310 

1290 

181 

4 

5 

1.1 

D 

305 

1170 

172 

3 

51 

0.8 

E 

219 

1085 

133 

5 

17 

1.1 

H 

189 

926 

119 

5 

42 

1.0 

P 

113 

988 

74 

5 

17 

0.6 

—  b 


Mean 


0.9 
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The  bedrock  on  Line  F  was  particularly  rough  and 
erratic  and  if  we  omit  this  low  value,  the  average  shear 
stress  at  the  bed  for  the  5  upper  lines  becomes  1.0  bars 
(1  bar  =  10°  dynes/cm  ).  Nye  (1952  b)  has  made  calcula¬ 
tions  for  sixteen  Alpine  glaciers  and  finds  values 
ranging  from  0.49  to  1.51  bars.  The  values  above  are 
in  the  middle  of  Nye’s  range  and  show  a  remarkably  small 
variation.  Glen  (1952)  has  used  uniaxial  compress! onal 
and  tensional  stresses  of  1.5  to  10.0  bars  in  the 
laboratory  to  produce  strain  rates  similar  to  those 
observed  on  glaciers. 

Now  that  we  have  an  average  value  for  the  shear 
stress  and  have  seen  that  it  does  not  depart  much  from 
a  value  of  1.0  bars,  we  can  make  use  of  relation  (12) 
to  determine  the  thickness  of  the  glacier.  The  value 
of  the  constant  hG  becomes  11.3  meters  or  37*1  feet. 
Formula  (12)  then  gives  a  depth  of  750  feet  for  D  line 
and  400  feet  for  E  line.  The  frictional  resistance  of 
the  sides  of  the  valley  has  been  neglected  in  this 
simple  expression  so  the  estimated  depth  is  too  small 
as  expected. 

In  point  of  fact,  formula  (12)  may  be  used  to 
calculate  the  thickness  of  an  ice  sheet  if  the  only 
information  available  is  the  surface  slope  (Nye,  1952  a). 
For  the  relation  to  be  satisfied,  the  sheet  of  ice  must 
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be  known  to  be  moving.  An  area  which  satisfies  this 
condition  is  the  Snow  Dome  on  the  Columbia  Icefields. 

It  has  an  average  slope  of  10°  to  11°  according  to  the 
Boundary  Survey  maps.  The  estimated  depth  of  ice  on 
this  feature  is  deduced  to  be  about  200  feet. 

Reference  to  Figure  13  discloses  that  the  glacier 
above  Line  B  has  a  steep  gradient  while  below  Line  B 
the  slope  is  very  small  and  even  becomes  negative.  If 
some  section  of  the  glacier  becomes  thicker  than  required 
by  equation  (12),  this  section  will  exert  a  pressure  on 
the  part  of  the  glacier  that  is  too  thin  and  relation 
(12)  will  be  re-established  by  flow.  According  to 
Orowan  (1949),  the  additional  longitudinal  force  required 
for  pushing  the  glacier  along  over  the  low  gradient 
part  may  be  so  large  that  it  exceeds  the  resistance  of 
the  ice  to  longitudinal  compression.  Longitudinal 
compression  of  the  glacier  takes  place  and  the  glacier 
thickens  until  the  cross  section  is  large  enough  to 
transmit  the  longitudinal  force  without  the  longitudinal 
yield  stress  of  the  Ice  being  exceeded. 

6 .  CONCLUSIONS 

The  thickness  of  the  Athabaska  glacier  has  been 

obtained  along  eight  transverse  lines  by  an  Investigation 
of  the  gravity  anomalies.  Previous  workers  have  obtained 
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the  depth  with  a  precision  that  varied  from  -10%  to 
+25%  or  in  some  cases  +40%.  A  higher  degree  of  accuracy 
was  attained  in  this  study  by  making  detailed  terrain 
corrections  and  by  treating  the  glacier  as  a  three 
dimensional  ore  body.  The  Bouguer  gravity  anomalies 

were  determined  within  0.5  milligals  or  better  and  the 

■» 

depth  values  are  believed  to  be  accurate  within  -10% 
to  +15%*  The  shape  of  the  bedrock  can  be  determined 
within  about  +  10%  if  thb  central  depth  is  known  from 
other  sources  and  the  edge  of  the  glacier  can  be 
determined  accurately. 

The  major  disadvantage  of  the  method  is  that  it 
requires  a  large  scale  contour  map  of  the  area  to 
reduce  the  terrain  corrections.  This  problem  is  rapidly 
being  overcome  for  most  mountainous  regions  by  the  use 
of  aerial  photogrammetry.  Making  the  terrain  corrections 
is  a  long  and  tedious  operation  but  this  can  be  over¬ 
come  by  use  of  Slegert’s  gravity  integrator  or  by  the 
use  of  a  digital  computer  with  a  large  storage  capacity. 

Previous  depth  determinations  were  limited  to  the 
smooth,  central,  portion  of  a  glacier  since  a  two 
dimensional  model  was  used  in  the  computations.  Three 
dimensional  computations  with  a  low  speed  computer  were 
made  in  this  Btudy  and  the  depth  near  the  terminus  of 
the  glacier  and  in  the  upper  elevated  section  was  obtained. 
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The  depth  of  the  Athabaska  glacier  varied  from  1070  feet 
on  line  B  to  160  feet  on  line  G. 

The  gravity  method  is  a  valuable  adjunct  to  a 
seismic  survey  since  a  detailed  set  of  measurements  can 
be  obtained  more  economically.  In  crevassed  areas 
where  seismic  reflections  are  difficult  to  obtain  and 
in  areas  where  multiple  reflected  events  are  recorded 
the  gravity  method  has  a  distinct  advantage. 

In  making  depth  computations ,  a  useful  first 
approximation  to  the  cross  section  of  the  valley  of  a 
glacier  is  a  parabola.  The  cross  section  of  the  valley 
through  which  the  Athabaska  glacier  flows  has  been 
found  to  approximate  a  parabola  along  several  of  the 
lines.  Svensson  has  noted  that  the  shape  of  some 
glacial  valleys  also  approximates  a  parabola.  It  is 
well  known  that  glaciers  alter  originally  stream-cut 
valleys  from  a  V  to  a  U-shaped  profile.  This  evidence 
would  seem  to  indicate  some  theoretical  significance  to 
a  parabolic  shape. 

A  basin  has  been  formed  In  the  center  of  the  glacier, 
probably  by  some  form  of  erosion  so  that  the  glacier  now 
flows  along  a  surface  that  is  flat  or  slopes  slightly 
upwards  toward  the  terminus.  Combined  gravity  and 
seismic  data  suggest  a  similar  up -valley  slope  on  the 
first  step  of  the  glacial  stairway.  If  the  glacier 


. 
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should  disappear  entirely,  the  chain  of  basins  would 
form  "paternoster"  lakes.  The  glacier  has  thickened 
to  its  maximum  extent  in  the  central  region  so  that  it 
now  has  a  surface  slope  which  allows  it  to  flow  along 
the  nearly  level  valley  floor. 

Prom  a  knowledge  of  the  depth,  shape  and  surface 
slope  of  the  glacier,  the  average  shear  stress  exerted 
by  the  bed  on  the  ice  was  found  to  be  1.0  bars.  *' 
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PART  II 

A  GRAVITY  INVESTIGATION  OF  THE 
ROCKY  MOUNTAIN  TRENCH  IN  SOUTHERN  BRITISH  COLUMBIA 
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1.  INTRODUCTION 

A  regional  gravity  investigation  of  the  Southern 
Cordillera  was  undertaken  by  Garland  and  Tanner  in  195^- • 
Gravity  measurements  obtained  in  the  Rocky  Mountain 
Trench  indicated  areas  with  a  considerable  thickness  of 
fill.  It  is  the  purpose  of  this  investigation  to 
examine  the  Trench  between  Cranbrook  and  the  International 
Boundary  with  a  more  detailed  set  of  gravity  measure¬ 
ments.  The  objective  was  a  more  accurate  delineation  of 
depressions  in  the  bedrock  which  are  filled  with  low 
density  fill  and  a  search  for  other  anomalies  associated 
with  faulting  and  intrusion.  An  interpretation  is  made 
of  a  gravity  lineation  discovered  along  the  Moyie  Lenia 


fault . 
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2.  GEOLOGY  OP  THE  AREA 

The  Rocky  Mountain  Trench  is  a  prominant  linear 
depression  separating  the  Rocky  Mountain  system  on  the 
east  from  the  Purcell,  Selkirk,  Monashee,  Cariboo, 

Omlneca  and  Cassiar  ranges  on  the  west.  (See  Figure  la.) 

It  extends  from  Flathead  Lake  in  Montana  to  the  British 
Columbia-Yukon  border.  The  tern  "Trench"  was  defined  by 
Daly  (1912)  as  a  "long,  narrow,  intermontane  depression 
occupied  by  two  or  more  streams  alternately  draining  the 
depression  in  opposite  directions." 

In  the  area  under  study  the  rocks  in  the  Purcell 
range  consist  of  over  30,000  feet  of  quartzites,  argillites 
and  limestones  called  the  Purcell  series.  (See  Table  la.) 
These  sediments  are  of  Pre Cambrian  Beltian  age.  Basalt 
flows  called  the  Moyie  intruslves  and  extrusions  of 
andesitic  lava  and  tuff  Interrupt  the  sedimentary 
sequence.  The  section  is  arched  up  to  form  the  north¬ 
ward  plunging  Purcell  geanticline.  Near  Cranbrook  the 
geanticline  is  cut  by  several  transverse  faults,  the 
most  important  of  which  is  the  Moyle  Lent a  fault. 

Kirkham  believes  the  Moyie-Lenia  is  an  overthrust 
fault  with  the  west  block  moving  upward  along  a  plane 
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TABLE  la.  Table  of  Formations 
(following  Leech  1958) 


ERA 

PERIOD  OR  EPOCH 

Thi ckness 

Formation  Lithology  (feet) 

CENOZOIC 

sand,  alluvium, 
till,  gravel 

MESOZOIC 

Lower  Cretaceous  ? 

monzonite, 

granodiorite 

Mississipplan 

Rundle  Gp.  limestone  2500 

Banff  silty  lime-  1300 

stone 

Exshaw  black  shale  57 

Devonian 

Palliser  massive  800 

GP*  limestone 

PALAEOZOIC 

Alexo  sandy  lime- n 

stone  ) 

\  1 335 

Fairholme  limestone  j 

Harrogate  black  lime¬ 
stone 

Bumais  gypsum, 

dolomite  600 

"Basal"  unit  sandy  dolomite 

+  +  +  +  + 

UNCONFORMITY  +  +  +  +  +  +  + 

Silurian  to 

Beaverfoot -Brisco,  Wonah,  Glenogle, 

Middle  Cambrian 

McKay,  Jubilee,  Elko,  and  Burton. 

Lower  Cambrian 

Eager  shale,  limestone,  ss. 

Cranbrook  quartzite 

+  +  *4-  + 

+  1*  +  +  4* 

UNCONFORMITY  +  +  +  +  +  +  + 

Purcell? 

Moyie  intrusion  meta-dlorite 

Upper  Purcell 

Dutch  Creek  doiomitic  3000 

argilites  &  ss. 

PRE CAMBRIAN 

Lower  Purcell 

Purcell  andesitic 

extrusions  lava  &  465 

tuff 

Siyeh  argillite 

Kitchener  doiomitic  > 7000 

argillites] 

Creston  argillites,  5000-6000 

quartzites 

Aldridge  argillites,  8000 

quartzites 

Fort  Steele  white 

quartzite  6000+ 
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never  less  than  45°  In  slope.  Rice  (1937)  states  that 
the  dip  Is  60°  west  at  Peavine  Creek  at  the  north  end 
of  Moyle  Lake.  The  Aldridge  formation  on  the  west  lies 
against  the  Siyeh  on  the  east  to  give  a  vertical  displace¬ 
ment  of  about  15,000  feet.  The  Moyie  Lenia  fault  has 
been  closely  positioned  by  Leech  in  the  Joseph  Creek 
area  south  of  Cranbrook  (written  communication,  see 
Figure  2a,  In  pocket ) . 

The  Purcell  series  occupies  the  entire  western 
edge  of  the  Rocky  Mountains  in  this  area.  The  broad 
structure  indicates  that  this  forms  the  east  flank  of 
the  Purcell  geanticline  but  locally  there  are  many 
small  folds  and  the  area  is  complicated  by  faulting. 

The  Dibble  Creek  fault  enters  the  Trench  transversely 
southeast  of  Fort  Steele.  Leech  (1959)  considers  that 
the  north  block  had  a  relative  upward  and  eastward 
movement  and  that  lubrication  was  provided  by  Devonian 
gypsum  (3urnais  formation)  in  the  footwall.  In  the 
Rocky  Mountain  ranges  northeast  of  the  Dibble  Creek 
fault,  Devonian  rocks  lie  on  Silurian,  Ordovician  and 
Cambrian  beds  but  south  of  the  fault,  the  Devonian  lies 
on  a  few  feet  of  Cambrian  or  directly  on  the  Kitchener- 
Siyeh  formation. 

The  floor  of  the  Trench  between  Bull  River  and  Sand 
Creek  contains  outcrops  of  Devonian  and  Mississlppian 
strata.  (See  Figure  4a).  Schofield  (1951)  found  these 
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to  rest  unconformably  on  the  Precambrian  Gateway  forma¬ 
tion  west  of  Wardner  at  Gold  Creek.  The  same  sequence 
of  Devonian -Mississippi an  rocks  reappears  in  the  Lizard 
range  south  of  the  Sand  Creek  Fault. 

Important  faults  have  been  traced  in  the  Trench 
between  Bull  River  and  Galloway  (Leech,  1958).  A  longi¬ 
tudinal  west  dipping  fault  separates  the  Purcell  series 
in  the  Rockies  to  the  northeast  from  the  Palaeozoics  in 
the  trench.  In  addition,  several  transverse  faults  cut 
the  Palaeozoic  strata.  Leech  (1959)  believes  that  the 
Palaeozoic  rocks  in  the  Trench  "belong  to  an  upper  plate 
dropped  along  the  west  dipping  fault  at  the  east  side 
of  the  trench  rather  than  that  they  belong  to  a  plate 
thrust  from  below."  He  also  finds  that  the  Precambrian 
rocks  on  the  trench  wall  are  themselves  a  thrust  plate 
and  partially  underlain  by  Palaeozoic  rocks. 

Several  small  outcrops  of  altered  monzonitic  or 
granodioritic  porphyry  indicate  a  small  stock  near  the 
mouth  of  the  Bull  River.  Sheep  Mountain,  west  of  the 
confluence  of  the  Elk  and  Wigwam  Rivers,  is  intruded  with 
igneous  rock  (written  communication.  Leech). 

The  trench  has  a  thick  veneer  of  glacial  till,  and 
indeed,  glacial  material  is  found  over  much  of  the  area 
except  on  the  highest  peaks.  The  rivers  have  cut  valleys 
several  hundred  feet  deep  through  the  material  in  the 
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Trench.  A  deposit  of  gravel,  sand  and  silt  known  as  the 
St.  Eugene  Silts  occurs  along  St.  Mary  River.  Schofield 
(1915)  regarded  these  as  interglacial  or  Pleistocene  in 
age  but  Hollock  and  Berry  consider  them  to  be  late 
Tertiary  from  a  study  of  fossil  plants.  If  this  is  so, 
they  somehow  escaped  erosion  by  the  ice.  Rice  (1937) 
points  out  several  facts  which  point  to  little  active 
erosion  during  the  period  of  glaciation.  These  include 
the  occurrence  of  V-shaped  tributary  valleys  and  placer 
gold  deposits  in  morainic  material. 

3.  THE  GRAVITY  MEASUREMENTS  AND  THEIR  REDUCTION 
3 • 1  The  Gravity  Survey 

Observations  were  carried  out  with  a  Worden  Gravity 
Meter  which  has  a  calibration  constant  of  0.2607  milligals 
per  scale  division.  About  300  gravity  stations  have 
been  established  in  the  Cranbrook-Elko  area.  One-third 
of  these  are  either  from  a  traverse  by  Garland  and  Tanner 
in  193^  or  by  Garland  in  1958-  T.  L.  Thompson,  a 
graduate  geology  student  from  Stanford,  who  was  studying 
the  structure  In  the  vicinity  of  Bull  River,  made 
observations  at  127  stations.  The  remainder  of  the 
observations  were  made  in  1959  and  i960  by  the  author. 

The  readings  were  taken  along  highways,  forestry 
trails,  railroads, or  at  stream  crossings  and  survey 
monuments.  The  location  of  the  stations  were  spotted  on 
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the  Columbia  River  Basin  maps  which  are  on  a  scale  of 
two  inches  to  a  mile  or  1:31*680.  The  contour  interval 
is  20  feet  up  to  an  elevation  of  3000  feet.  Maps  on  a 
scale  of  1:40,000  or  1:50,000  with  100  foot  contour 
interval  were  used  for  the  remainder  of  the  area.  All 
gravity  measurements  were  made  with  respect  to  the 
network  of  stations  established  by  the  Dominion  Observa¬ 
tory  (Garland  and  Tanner,  1954).  Many  ties  were  made 
between  the  stations  obtained  at  different  times  and 
the  results  check  within  0.2  milligals.  A  direct 
comparison  can  be  made  at  stations  6,  10,  13,  23,  33, 
lol,  166,  134,  204,  203  and  211. 

Gravity  and  aneroid  barometer  readings  were  made  in 
loops  from  a  base  station  and  repeated  at  least  once 
every  two  hours.  A  linear  drift  rate  was  assumed  for 
both  meter  and  barometer  and  the  readings  at  new  stations 
were  corrected  for  this  drift.  The  barometer  readings 
were  used  whenever  station  elevations  were  not  available 
from  bench  marks  and  spot  highv;ay  elevations.  The 
elevation  of  most  of  the  stations  obtained  in  this  manner 
are  believed  accurate  within  5  feet  but  some  may  be  as 
much  as  10  feet  off.  However,  one  line  of  stations 
(Numbers  47  to  58  inclusive)  is  in  douot  by  +  35  feet. 

The  Bouguer  values  on  this  line  may  be  in  error  by  +  2.1 
milligals.  The  station  locations  are  believed  accurate 
within  about  500  feet. 
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3 . 2  The  Reduction  of  the  Gravity  Measurements 

The  principal  facts  for  the  gravity  stations  are 
set  out  in  Table  Ila.  This  table  gives  the  station, 
field  number,  longitude,  latitude,  elevation,  observed 
gravity  and  the  corrections  necessary  to  obtain  the 
Bouguer  anomaly.  The  theoretical  gravity  at  sea  level,  y , 
was  obtained  from  the  1930  International  Gravity  Formula 
(Nettleton,  19^0).  The  free  air  and  Bouguer  correction 
amounts  to  0.060  milligals  per  foot  above  sea  level.  A 
density  of  2.6?  gms.  per  cc.  was  used  in  correcting  for 
the  attraction  of  material  between  the  station  and  sea 
level.  This  is  consistent  with  the  practice  followed 
by  the  Dominion  Observatory  and  is  close  to  the  density 
of  Devonian  and  Carboniferous  rocks  in  this  area.  If 
Ii  Is  the  elevation  of  the  station  in  feet  above  sea 
level,  the  simple  Bouguer  anomaly  is  given  by  the  ex¬ 
pression  : 

6  -  sob8erved  -  7  +  °-o6oH 

Time  did  not  permit  the  application  of  terrain 
corrections  to  the  measurements.  This  type  of  correction 
Is,  of  course,  desirable  for  a  full  interpretation  of 
any  mountainous  region.  However,  Garland  and  Tanner 
(195^0  have  made  terrain  corrections  to  53  of  the 
stations  used  In  this  map  area.  An  examination  of  these 
corrections  shows  that  they  are  small  and  uniform  over  the 
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area  of  the  trench.  An  understanding  of  this  uniformity 
can  be  gained  fr^m  a  study  of  the  physiography  of  the 
area  (Figure  la).  Nearly  all  the  measurements  were 

taken  in  the  wide,  fairly  level,  depression  making  up  the 
Trench.  Nearly  all  the  stations  have  an  elevation 
within  the  range  between  2500  and  3500  feet.  The  Purcell 
range  has  an  average  elevation  of  5000  feet  with 
occasional  peaks  rising  to  7000  feet  above  sea  level. 

The  mountains  merge  gradually  into  the  trench  but 
isolated  peaks  such  as  Mt.  Baker  do  occur  on  the  edge. 

The  Rocky  Mountains  rise  in  an  abrupt  wall  6000  to  8000 
feet  above  sea  level.  South  of  Cranbrook  the  Trench  is 
generally  8  to  14  miles  wide.  The  gravitational  attrac¬ 
tion  of  the  Rocky  Mountain  mass  is  very  nearly  the  same 
on  most  of  the  stations  in  this  wide  depression.  This  is 
best  illustrated  in  Table  Ilia  summarizing  the  terrain 
correction  made  by  Garland  and  Tanner. 

TABLE  Ilia.  Terrain  Corrections  in  Cranbrook- 
Elko  Area 


Region 

Terrain  Correction 
(milligals ) 

Average 
Correction 
(mi lligals ) 

Number 

of 

Stations 

Rocky  Mountain 

0.5  -  1-5 

26 

Trench 

1.6  -  2-5 

2.6  -  3-5 

1.4 

9 

6 

(41) 

Purcell  Mountains 

1.5  -  3.8 

2.7 

7 

Rocky  Mountains 

3-4  -  9.3 

_ 

6.  o 

5 

75 

The  conclusion  is  reached  that  the  effect  of  terrain 
is  small  and  uniform  along  the  Trench  and  the  valley  of 
the  St.  Mary  River.  It  is,  therefore,  possible  to 
treat  the  terrain  effects  as  part  of  the  regional  gradi¬ 
ent  in  the  interpretation  of  the  gravity  anomalies. 

3. 3  Density  of  the  Rocks 

Garland  and  Tanner  found  a  density  difference  of 
0.7  gms.  per  cc.  between  the  unconsolidated  material 
along  the  Columbia  River  and  bedrock.  This  high  value 
can  only  apply  to  surficial  deposits  and  would  be 
expected  to  decrease  as  the  fill  was  buried  deeply.  If 
the  St.  Eugene  silts  are  Miocene  as  concluded  by  Berry, 
there  is  a  good  chance  of  other  Tertiary  deposits  in 
the  trench.  The  density  difference  between  the  Tertiary 
deposits  and  Palaeozoic  or  Pre Cambrian  bedrock  might  be 
expected  to  be  as  low  as  0.3  or  0.4  gm.  per  cc.  Since 
there  are  no  measurements  of  density  within  the  trench, 
a  value  of  0.5  gms.  per  cc.  will  be  used  in  the  inter¬ 
pretation  of  anomalies  within  the  Trench.  The  calculated 
depth  to  bedrock  is  inversely  proportional  to  the  density 
difference  so  a  reasonable  estimate  of  the  accuracy  can 
be  made. 

Deposits  of  gypsum  occur  in  the  Devonian  Burnais 
formation.  Garland  and  Tanner  have  found  their  density 
to  be  2.26  gm.  per  cc.  These  deposits  occur  along  the 
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Dibble  Creek  fault  and  along  Bull  River.  These  features 
would  be  expected  to  show  up  as  narrow  bands  of  negative 
gravity  anomalies. 

Small  Intrusive  bodies  with  compositions  ranging 
from  granodiorite  to  syenite  occur  north  of  St.  Mary 
.River,  at  the  mouth  of  Bull  River  and  in  Hughes  Range. 

The  density  of  this  material  is  not  known  but  there  is 
a  small  positive  anomaly  over  the  stock  of  altered 
monzonitic  or  granodioritic  porphyry  near  the  mouth  of 
the  Bull  River. 

Lastly,  the  region  of  the  Moyie  Lenia  fault  will  be 
considered.  The  nature  of  the  fault  is  not  too  important 
to  gravity  interpretation  as  the  formations  on  both 
sides  are  similar  and  would  have  a  similar  density. 
However,  the  Moyie  intrusives  formed  dense  sills  in  the 
Aldridge  formation  and  several  lava  flows  occur  near  the 
top  of  the  Siyeh  formation.  These  might  be  expected  to 
provide  sufficient  density  contrast  to  delineate  the 
fault  zone.  The  Purcell  lavas  vary  in  composition  from 
an  andesite  to  a  diabase.  On  Mt.  Baker  they  vary  in 
thickness  from  50  to  300  feet  and  each  flow  is  separated 
by  50  to  400  feet  of  sediment  (Schofield,  1915).  Garland 
(1954)  finds  the  mean  density  of  the  Lower  Purcell  series 
to  be  2.74  gm.  per  cc.  A  sample  of  amygdaloidal  lava 
near  Skookumchuck  had  a  density  of  2.84.  Reich  found 
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diabase  extrusives  to  have  densities  varying  from  2.73 
to  3.12.  A  density  difference  of  0.10  gras,  per  cc. 
between  the  sediments  and  the  extrusion  is  probably  close 
to  the  true  value. 

The  Moyie  intrusives  vary  in  thickness  from  2  to 
over  2000  feet.  Schofield  refers  to  these  as  the  Purcell 
sills  and  Rice  calls  them  the  Purcell  intrusives.  Their 
composition  varies  from  a  gabbro  to  a  quartz  diorite 
to  granite.  The  densities  of  two  typical  sections  was 
determined  by  Schofield  (1915)  and  is  reproduced  in 
Table  IVa.  By  weighting  each  density  measurement 
according  to  the  thickness  the  mean  density  of  the 
gabbro  is  found  to  be  2.99  gms .  per  cc.  and  the  mean 
density  of  the  granite  is  2.75*  The  mean  weighted 
density  of  the  3820  foot  column  at  Kingsgate  is  2.39  gms. 
per  cc.  and  the  mean  weighted  density  of  the  5226  foot 
column  at  St.  Mary  Lake  is  2.93  gm.  per  cc.  It  would 
seem  that  there  is  a  density  excess  of  about  0.20  gms. 
per  cc.  between  a  purely  sedimentary  sequence  and  one 
intruded  with  Moyie  sills. 

3 . 4  Accuracy  of  the  Data 

The  simple  Bouguer  anomalies  are  believed  to  be 
accurate  within  0.2  milligals  for  all  stations  located 
at  benchmarks  or  surveyed  points  along  highways  and 
railroads.  The  effect  of  location  errors  are  negligible 
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Table  IVa.  Densities  of  Two  Columnar  Sections  In  the  Purcell 

Range 

(after  Schofield,  1915) 


Formation 

Rock  Type 

Thi ckness 

Average  Density 

9 

(feet) 

(gm/cc. ) 

1.  Moyle  Sill  on  Mountain  west  - 

of  Kinsgate 

Aldridge 

quartzite 

200 

2.76 

Moyle  Sill 

A 

gabbro 

2  6 

2.96 

granite 

80 

2.76 

gabbro 

29 

2.97 

Aldridge 

quartzite 

100 

2.766 

Movie  Sill 

B 

gabbro 

30 

2.93 

Aldridge 

quartzite 

670 

2.76 

Moyle  Sill 

C 

granite 

200 

2.74 

intermediate 

87 

2.76 

gabbro 

623 

3-00 

Aldridge 

quartzite 

75 

2.76 

Moyie  Sill 

D 

gabbro 

1500 

2.99 

Aldridge 

quartzite 

200 

2.76 

Total  Thickness 

3820 

2.  St 

; . 

Mary  Sills  on  Mountain  10  miles  west 

of  Marysville 

Aldridge 

quartzite 

Sill  A 

granite 

70 

2.76 

gabbro  &  intermediate 

70 

2.89 

Aldridge 

quartzite,  argillaceous 

400 

Sill  B 

gabbro  (hornblende) 

985 

2.99 

Aldridge 

quartzite,  argillaceous 

200 

Sill  C 

gabbro,  hornblende 

123 

2.99 

Aldridge 

quartzite 

300+ 

Sill  D 

gabbro 

565 

2.99 

Aldridge 

quartzite 

350+ 

Sill  E 

gabbro 

2l63 

2.99 

Total  thickness 

5226 
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as  a  mistake  of  1000  feet  in  latitude  must  be  made  to 
introduce  0.1  milligals  error  in  the  gravity  value.  The 
precision  of  stations  for  which  the  elevations  were 
obtained  by  aneroid  barometer  is  +  0.6  milligals  except 
for  stations  47  to  58  which  are  in  doubt  by  +  2.1  milligals. 

4.  INTERPRETATION  OF  GRAVITY  ANOMALIES 

4 . 1  Anomalies  Due  to  Overburden 

Gravity  methods  are  very  effective  in  determining 
the  thickness  of  unconsolidated  material  filling  a  basin 
around  which  bedrock  outcrops.  While  the  exact  shape  and 
depth  to  bedrock  is  indeterminate,  the  approximate 
result  is  often  enough  to  determine  the  origin  of  the 
basin,  be  it  glacial  erosion  or  normal  faulting  or  some 
other  mechanism.  The  Interpretation  can  proceed  in  a 
straightforward  manner  because  one  face  of  the  body 
producing  the  anomaly,  that  is  the  surface,  is  completely 
determined.  The  density  contrast  between  the  fill  and 
bedrock  is  so  large  that  small  irregularities  in  the 
densities  do  not  affect  the  overall  conclusion.  The 
results  of  two  investigations  of  similar  anomalies  in 
Nevada  and  California  have  recently  been  published 
(Thompson  and  Sandberg,  1958;  Pakiser  et_  al,  i960). 

In  many  geologic  situations,  the  horizontal  dimen¬ 
sions  of  the  density  discontinuities  are  several  times 
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the  vertical  dimensions.  In  such  a  case,  the  thickness 
of  the  body  can  be  approximated  by  the  Bouguer  expression 
for  an  infinite  plate. 

A  q  -  2  77-  k cr  H  (i) 

K  =  the  gravitational  constant 
G~  =  density 

=  thickness  of  the  plate 

If  the  density  difference  between  bedrock  and  fill  is 
0.5  gms.  per  cc.  the  thickness  of  the  deposit  in  feet  is 
approximately : 


H  =  157  Aq  (2) 

Because  the  body  has  finite  horizontal  dimensions,  the 
depth  will  be  slightly  underestimated. 

If  the  gravity  anomaly  is  roughly  circular  the  body 
can  be  approximated  by  one  or  more  cylinders.  The 
gravity  anomaly  at  the  center  of  the  cylinder  can  be 
found  exactly  (Heiland,  1940,  p.  147).  The  approximate 
size  and  thickness  of  the  cylinders  can  be  obtained  from 
an  application  of  the  Bouguer  expression  (2).  The 
gravity  anomaly  of  a  surface  distribution  of  mass  is 
proportional  to  the  solid  angle  subtended  by  the  surface 
distribution  at  the  station  in  question.  Nettleton 
(1942)  has  tabulated  the  solid  angles  subtended  by 
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circular  laminar  discs.  The  mass  of  the  cylinder  is 
then  considered  to  be  condensed  upon  the  median  plane 
and  the  gravity  anomaly  obtained  from  the  tabulated 
solid  angles.  It  is  then  possible  to  calculate  the 
approximate  anomaly  at  stations  off  the  origin. 

An  examination  of  the  Bouguer  map  reveals  three 
large  negative  anomalies  (Figure  2a,  in  pocket).  The 
amplitude  of  the  anomalies  decrease  to  the  northwest. 

The  largest  and  southern -most  anomaly  is  the  Waldo 
gravity  minimum  (Figure  3a).  North  of  the  Waldo 
minimum  is  a  small  positive  anomaly  centered  at  Elko. 

It  is  underlain,  in  part,  by  the  Dutch  Creek  formation 
around  Elko,  and  in  part,  by  intrusives  outcropping  at 
Sheep  Mountain.  This  high  represents  the  higher  density 
of  the  Purcell  series  (2.7^  gras,  per  cc.)  as  compared 
to  the  density  used  in  making  the  Bouguer  correction 
(2.67)»  From  longitudinal  and  transverse  cross  sections 
of  the  Bouguer  anomalies,  the  regional  gradient  can  be 
estimated.  This  regional  gradient  Is  due,  in  part,  to 
terrain  effects,  and  in  part,  to  deep  seated  density 
variations  (isostatic  compensation)  which  are  not  of 
interest  here.  The  residual  anomaly,  relative  to  out¬ 
crops  of  bedrock,  amounts  to  -19  milligals  on  this 
southern-most  minimum.  If  the  density  difference  between 
fill  and  bedrock  is  -0.5  gms.  per  cc.,  then  the  thickness 
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of  fill  is  3500  feet.  The  density  difference  is  thought 
to  have  a  precision  of  +0.1  gms.  per  cc.  The  calculated 
depth  may  be  as  great  as  4400  feet  or  as  small  as 
2900  feet . 

The  fourteen  milllgal  gravity  minimum  at  Jaffray 
(Figure  4a)  is  in  the  shape  of  a  square,  three  miles  to 
a  side.  The  depth  of  fill  is  calculated  to  be  2300 
feet.  From  scattered  outcrops  to  the  north  and  east 
of  this  anomaly  Leech  (1958)  maps  a  west  dipping  fault 
on  the  east  wall  of  the  Trench.  There  is  a  large 
gravity  gradient  which  confirms  the  presence  of  this 
fault.  The  area  north  of  this  fault  is  particularly 
interesting  as  it  consists  of  a  Palaeozoic  sequence 
that  rises  2200  feet  above  the  level  of  the  trench  floor. 
It  is  entirely  broken  up  by  faulting  and  is  evidently 
a  small  horst  adjacent  to  the  graben  making  up  the 
Jaffray  gravity  minimum. 

The  Fort  Steele  gravity  minimum  (Figure  4a)  has 
eight  milligals  of  closure.  It  is  rectangular  in 
shape  with  dimensions  three  by  four  miles.  The  calculated 
depth  of  fill  is  1400  feet.  The  northern  edge  is 
complicated  by  the  presence  of  the  Boulder  Creek  fault. 
This  fault  appears  to  be  a  thrust  with  the  Aldridge 
formation  on  the  north  and  the  Creston  and  Kitchener 
formations  on  the  south  but  there  may  be  a  strike- 
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Fif^ure  4a.  Gravity  Profile  and  Geololc  Goction  across  the  Roclcy 
Mountain  Trench  at  Fort  Gteole  and  Jaffray.  The  locations  of  the 
profiles  ore  shown  in  figure  2a.  The  horizontal  scale  is  twice 
the  vortical  scale. 
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slip  component  to  the  motion.  The  dip  of  the  fault  is 
very  nearly  vertical  or  steeply  northward  and  it  brings 
together  blocks  with  a  complicated  history  (Leech,  1958). 

This  fault  may  extend  to  the  St.  Mary  fault  in  the 
Purcell  Mountains  but  more  gravity  measurements  are 
necessary  to  establish  this  correlation. 

North  of  Fort  Steele  the  bedrock  must  be  beneath 
a  thin  cover  of  till.  There  is  a  suggestion  of  a 
fourth  minimum  four  miles  north  of  Fort  Steele.  The 
Dibble  Creek  fault,  which  has  gypsum  in  the  foot wall,  may 
extend  into  the  trench  and  widen  the  southern  end  of 
the  Fort  Steele  minimum. 

4.2  Anomalies  Related  to  Faults  and  Intrusives 

Various  small  anomalies  are  apparent  on  the  Bouguer 
map  which  appear  related  to  structural  features.  The 
density  of  the  data  is  not  sufficient  to  make  quantitative 
calculations  but  a  qualitative  discussion  may  stimulate 
further  research.  Any  gravity  study  which  is  to  be 
successful  depends  critically  on  a  prior  detailed 
geologic  study.  The  density  contrast  between  till  or 
alluvium  and  bedrock  is  so  great  that  gravity  anomalies 
due  to  the  recent  deposits  may  entirely  obscure  more 
fundamental  trends. 

The  stock  of  monzonitlc  or  granodloritic  porphyry  at  Bull  River 
appears  to  cause  a  one  mllligal  gravity  high.  Intrusives 
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on  Sheep  Mountain  south  of  Elko  appear  to  cause  an  anomaly 
of  similar  magnitude.  The  composition  of  the  lntrusives 
on  Sheep  Mountain  is  not  known.  If  this  correlation  is 
correct  the  bulk  of  the  intrusive  must  have  a  composition 
which  corresponds  to  the  "grey  granodioritic"  phase  of 
Rice  (19^7)*  This  phase  has  very  little  potash  feldspar1 
but  about  14  per  cent  of  ferromagnesian  minerals.  The 
density  excess  over  the  surrounding  Palaeozoic  rocks 
can  scarcely  be  much  more  than  0.05  gms.  per  cc. 

The  low  trend  along  Bull  River  to  Lime  Creek  is 
probably  due  to  gypsum  which  has  a  density  of  2.26  gms. 
per  cc.  Further  north  the  anomaly  becomes  positive  so 
the  deposit  of  gypsum  may  follow  Lime  Creek  as  inferred 
by  Leech  (1958)  from  the  presence  of  sinkholes.  The 
gravity  anomaly  between  stations  31  and  8l  may  be  an 
extension  of  this  gypsiferous  formation  along  the  west 
(lipping  fault  mentioned  in  connection  with  the  Jaffray 
minimum. 
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4. 3  The  Moyle  Lenla  Fault 

A  linear  gravity  minimum  trends  along  the  Moyie 
Lenia  fault  south  of  Cranbrook  to  Moyle  Lake  (Figure  2a, 
in  pocket).  This  minimum  does  not  appear  to  be  related 
to  low  density  fill  because  several  low  values  are  near 
bedrock.  The  Purcell  sediments  on  both  sides  of  the 
fault  have  nearly  the  same  density  so  they  could  not 
materially  affect  the  gravitational  field.  It  is  pro¬ 
posed  that  the  Moyie  sills  are  responsible  for  the 
minimum  trend.  'The  sills  of  predominantly  basic  material 
have  intruded  the  Aldridge  formation.  The  density 
relations  have  already  been  reviewed  in  section  3»3- 

Figure  5a  illustrates  the  gravity  anomaly  across  the 
Moyie  Lenia  fault.  3elow  it  is  a  cross  section  showing 
two  solutions  which  produce  a  similar  gravitational 
field.  The  theoretical  gravity  produced  by 'the  Moyie 
sills  and  Purcell  extrusions  has  been  calculated  with  the 
aid  of  a  two  dimensional  graticule  chart  (Nettleton,  194o)‘. 
More  elaborate  computations  are  scarcely  warranted  until 
more  complete  geologic  and  geophysical  data  is  available. 
The  Moyie  Lenia  fault  presents  one  of  those  unusual 
situations  in  which  a  density  excess  produces  a  negative 
gravity  anomaly.  A  closer  fit  between  the  observed  and 
theoretical  gravity  could  be  made  by  altering  the  thickness 
and  dip  of  the  Moyie  intrusives  and  by  adding  a  small 
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Figure  5?.  Gravity  Profile  and  Geologic  Section  across  the 
Moyie  Lenia  Fault  in  the  iurcell  Mountains.  The  location  of 
the  profile  is  s.iown  in  f  ure  2a.  The  horizontal  and  vertical 
scales  are  e  ,ual.  Two  possible  inter  rotations  are  illustrated 
together  with  the  resultant  gravity  anomaly. 
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amount  of  fill.  Intense  shearing  of  the  rock  along  the 
fault  plane  may  have  decreased  the  density.  On  Peavine 
Creek  the  Aldridge  sediments  and  sills  are  converted 
from  massive  rocks  to  closely  cleaved,  fissile  phyllitic 
rocks  in  many  places.  (Leech,  written  communication). 

Much  of  the  area  is  blanketed  with  a  heavy  cover  of 
timber  and  outcrops  are  obscurred  by  recent  deposits  so 
the  thickness  of  the  section  containing  the  sills  is 
not  known.  In  the  theoretical  cross  section  it  is 
shown  to  increase  in  thickness  from  2000  to  5000  feet 
towards  the  southeast.  The  interpretation  is  not 
changed  significantly  if  the  single  sill  shown  in  the 
diagram  is  actually  composed  of  numerous  basaltic 
flows  separated  by  sediment. 

The  Purcell  extrusives  on  top  of  the  Siyeh  formation 
contribute  a  small  positive  anomaly  on  the  southeast 
block.  These  extensive  lavas,  mostly  dipping  20°  to 
30°  northeasterly,  may  total  about  1000  feet  in  thickness. 
The  section  is  repeated  by  faults,  most  of  which  appear 
to  be  steep  and  to  strike  northwesterly.  As  a  result 
the  plan  widths  of  lava  are  anomalously  great  (written 
comrauni cation.  Leech).  It  is  difficult  to  calculate  the 
net  effect  of  these  lavas  without  more  exact  geologic 
data  but  they  would  be  expected  to  produce  a  positive 
anomaly  of  one  or  two  milllgals. 
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The  normal  fault  solution  produces  a  better  fit  but 
this  solution  is  not  unique  and  a  slight  redistribution 
of  the  mass  in  the  overthrust  case  would  also  produce 
fair  agreement.  Kirkham  (1930 )  reviewed  all  the  data 
on  this  fault  and  concluded  that  it  was  an  overthrust 
with  a  steep  fault  plane  that  dips  west.  The  fault 
plane  forms  an  arc,  concave  to  the  east,  and  runs  for 
118  miles  through  British  Columbia,  Idaho  and  Montana. 

Daly  (1910)  believed  the  fault  to  be  normal  at  the 
International  Boundary.  Schofield  (1915)  states: 

The  second  great  movement  which  affected  the 
Purcell  Range  was  normal  faulting,  a  northeast  - 
southwest  system,  including  the  Marysville  and 
Moyie  faults  and  several  small  faults  of  lesser 
importance;  and  a  northwest -southeast  system, 
the  chief  one  being  the  Cranbrook  fault. 

However,  his  structure  diagram  (AB,  Map  147A)  shows  the 
Moyie  fault  to  be  overthrust.  The  problem  has  not  been 
resolved  in  this  study  but  it  appears  that  a  detailed 
gravitational  and  magnetic  survey  along  the  entire 
length  of  the  break  may  ultimately  provide  enough  data 
to  decide  if  the  fault  is  normal,  reverse,  or  possibly 
strike-slip. 

The  relative  low  along  the  Moyle  Lenia  fault  diminishes 
as  the  trench  is  approached.  This  is  readily  under¬ 
standable  as  the  east  limb  of  the  Purcell  geanticline 
dips  steeply  to  the  east  along  the  trench.  The  Siyeh 
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formation  outcrops  near  Rampart  and  this  indicates  that 
the  Moyie  intrusives  in  the  northwest  fault  block  have 
been  buried  too  deeply  to  cause  a  detectable  gravity 
anomaly.  It  is  doubtful  if  a  direct  connection  can  be 
shown  between  the  Moyie  Lenia  and  the  Dibble  Creek  fault 
by  gravity  methods. 


5.  CONCLUSIONS 

Gravity  measurements  indicate  the  presence  of  three 
deep  basins,  now  filled  with  unconsolidated  material, 
along  the  Trench  between  Fort  Steele  and  the  International 
Boundary.  From  north  to  south  the  depth  of  the  basins 
are  1400  feet  (Fort  Steele  minimum),  2300  feet  (Jaffray 
minimum),  and  3500  feet  (Waldo  minimum)  with  the  limits 
being  about  -11%  to  +25$.  The  basins  are  thought  to 
represent  down-faulted  blocks  along  longitudinal  and 
transverse  faults  along  the  trench.  Between  the  Jaffray 
and  Fort  Steele  minimum  there  is  2200  feet  of  topo¬ 
graphic  relief  with  Palaeozoic  outcrops  which  indicate 
that  this  is  a  small  horst.  Leeck(l959)  has  observed  that 
"South  of  latitude  49°  30',  the  Trench  seems  due  funda¬ 
mentally  to  block  faulting  and  lacks  important  strike- 
slip  displacements."  The  gravity  data  confirms  this 
observation  and  provides  evidence  that  the  amplitude  of 
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block  faulting  decreases  as  one  proceeds  northward  from 
the  International  Boundary  towards  Fort  Steele. 

This  investigation  has  indicated  an  economic 
significance  to  the  use  of  gravity  in  this  area.  Stocks 
of  raonzonitic  or  granodioritic  porphyry  appear  to  cause 
one-milligal  gravity  highs.  There  Is,  therefore,  a 
possibility  of  mapping  extensions  to  the  outcrops  with 
a  very  accurate  and  detailed  gravity  survey.  Many  of 
the  sulfide  deposits  in  this  area  occur  in  Aldridge 
strata.  Because  of  the  presence  of  dense  Moyie  sills  the 
appearance  of  the  Aldridge  formation  is  accompanied  by 
a  positive  gravity  anomaly.  This  is  Illustrated  between 
Marysville  and  Kimberley,  north  of  Fort  Steele  and  near 
Cranbrook.  Of  more  direct  application  Is  the  use  of 
gravity  to  determine  the  trend  and  Importance  of  gypsum 
deposits  which  are  used  in  the  cement  industry. 

A  linear  gravity  minimum  has  been  found  along  the 
Moyie  Lenia  fault.  This  trend  is  really  composed  of 
two  high  anomalies,  one  on  each  side  of  the  fault,  which 
are  caused  by  dense  sills  of  basaltic  material.  This 

f 

fault  can  be  treated  as  a  steeply  dipping  overthrust 
in  which  the  northwest  block  moved  upward  relative  to 
the  southeast  block  or  as  a  normal  fault  with  the  north¬ 
west  block  moving  upwards  relative  to  the  southeast 
block.  Of  the  two  Interpretations  presented  the  normal 
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fault  solution  fits  the  observed  gravity  data  slightly 


better.  It  does  not  seem  possible  to  tie  the  Dibble 

f 

Creek  fault  to  the  Moyie  Lenla  fault  on  gravity  data 
alone  but  the  Dibble  Creek  fault  appears  to  extend  into 
the  trench  to  within  three  miles  of  Rampart.  Similarly 
the  Boulder  Creek  fault  would  seem  to  extend  into  the 
Trench  past  Fort  Steele.  More  observations  are  necessary 

to  check  its  connection  to  the  St.  Mary  fau^t  system  as 

£ 

postulated  by  Rice  (1937)* 
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0s  O'  O' 
XXX 


O  <"Acm 
A-X  A- 
O'  O'  O' 
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X  -d  LA 
O'X  _d 
O'ccco 
XXX 


CM  O'  <A 

(Adx 

XXX 

XXX 


x  O'-d 

A-X  X 
XXX 
XXX 


© 

rd 

d 

A-  o  CM 

LAX  O 

rH  (ALA 

X  X  rH 

(A-dX 

C  O'  o 

o  o  o 

O  O  rH 

i — 1  rH  rH 

rH  O  rH 

i — 1  rH  i — 1 

i — 1  rH  AJ 

•p 

CM  <ACM 

CM  CM  CM 

CM  CM  CM 

CM  CM  AJ 

CVJ  CM  CM 

CM  CM  CM 

•H 

rH  rH  rH 

i — 1  i — 1  i — 1 

rH  rH  rH 

rH  rH  r— 1 

rH  rH  rH 

i — 1  < — 1  i — 1 

© 

>d 

CM  CM  CM 

CM  CM  CM 

CM  CM  CM 

CM  CM  CM 

CM  CM  CM 

CM  CM  CM 

o 

LA  LA  LA 

LA  LA  LA 

LA  LA  LA 

LA  LA  LA 

LA  LA  LA 

LA  LA  LA 

© 

d 

O  CM  CM 

<A_d-A 

X  A-  A- 

X  o  O 

rH  CM  CM 

<A_d-_d 

- 

la  _d_d 

-d'-d-d 

-d_d-d 

-d_d^t 

_d_rt-d 

-d  4  4 

rH  i — 1  i — 1 

i — 1  i — 1  r-H 

rH  rH  i— 1 

rH  rH  rH 

r-i  r-i  r~i 

i — 1  rH  i — 1 

bJO 

C 

A-  A- A- 

a-  r —  a~ 

A- A- A- 

A-  A-  A- 

A-  C —  A- 

A- A- A- 

O 

i — i  rH  i — 1 

i — (  i — 1  rH 

i — 1  rH  < — 1 

(Hi — 1  rH 

i — 1  i — 1  rH 

i — 1  rH  i — | 

id 

o 

H  H  H 

rH  i — 1  rH 

i — 1  I — 1  i — 1 

r-^  f-i 

i — 1  rH  I — 1 

rH  rH  rH 

e 

o 

O 

•rl 

1 

-P 

O'  rH  CM 

CA_dLA 

X  A-X 

O'  O  H 

CM  (Ad 

LAX  A- 

© 

P 

CO 

WdK 

K  ffi  W 

W  tH  tii 

Ph 
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4 

o 

LAX  A- 

X  O'  o 

rH  CM  (A 
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O  rH  CM 

S5 

LA  LA  LA 

LA  LAX 
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W 

H 

O 


O 

a 

CO 

m 
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Bouguer 

Anomaly 

milligals 

-166.9 

-165.6 

-182.6 

-178.0 

-176.1 

-178.6 

-4  CM  LA 
•  *  • 

CO  CO  CO 

r— r-r- 

HHH 

i  i  i 

-178.6 

-178.7 

-191.1 

CO  roro 
•  •  • 

rHLAS3 

O'  O'  C' 

rH  r — 1  rH 

1  1  1 

-196.7 

-197.3 

-197.6 

Terrain 

Dorrection 

(A-X) 

milligals 

S3  O'  C\J 

•  •  • 

r— un  c\j 
*  •  • 

co  co  co 

O'OCO 

•  *  • 

r-co  r- 

OV\cr\ 

9  •  • 

coco  o 

rH 

S3  COLA 

®  •  • 

CO  CO  O' 

©  CO  O' 

•  •  • 

O'  O' CO 

00 

rH 

©  (>> 

© 

W\00' 

rH  CO  r- 

coco  co 

s)  _4r- 

0s  co  r-» 

Oco-4 

rH  ©  rH 

bO 

LALAf— 

r—tAr- 

00  rH  CA 

LA  CO  CO 

r-His 

cA4"iA 

a  3  cc 

•H 

•  •  • 

•  •  • 

9  4  9 

•  •  • 

•  •  9 

•  4  • 

E  to  g 

< - 1 

<A  CM  O' 

vO  S3  S3 

S3  S3  S3 

vOt —  rH 

O  IA1A 

S3  S3  S3 

■H  3  O 

f - 1 

O'  O' 33 

CO  CO  CO 

coco  GO 

COCO  O 

OOO 

o  o  o 

woe 

•H 

i — 1  i — 1  rH 

rH  f— 1  i — 1 

t — 1  i — 1  rH 

rH  i — 1  OJ 

co  co  co 

co  CO  CO 

CQ  < 

E 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

i  ^ 

CQ 

•H  >>> 

rH 

I 

© 

CO  rH  S3 

CO  OOO 

CABALA 

CO  44 

OJ  rH  rH 

s 

,b£) 

CO  LA  ca 

-Cj-c^rH 

LA  S3  O 

roroo 

0-lA  rH 

O  CM  S3 

®  e 

«rl 

©  o 

rH 

O  H  CO 

r-co  go 

co  co  co 

S3  S3_4 

4ocr 

co  r-sD 

a 

rH 

-4-4-4 

rA^ArA 

cArA  ca 

CO  cOvO 

S3  LA -4 

fc,  < 

•H 

i 

S 

m 

OOIA 

cac—o 

CO  co  O 

OJ  O  co 

0s 'LA  O 

r-CO  rH 

©  t>5 

rH 

or- 

O  co  -o 

COLOCO 

LA  CO  rH 

S3  O  03 

COCO  <A 

I  >  -p 

© 

bD 

sO  CO  c\] 

CO  o  o 

rH  i — 1  OJ 

_4co  co 

CO  S3  S) 

c— go  O' 

©  > 

•H 

C'O'O' 

O'  O'  O' 

O''  0s  O 

O  rH  rH 

rH  rH  i— 1 

to  a) 

rH 

vOvO^O 

S3  O  S3 

©©© 

S3  S3  S3 

©  ©  © 

S3  S3  S3 

&  t* 

rH 

O  c3 

oH 

o  o  o 

OOO 

OOO 

OOO 

OOO 

o  o  o 

E 

23  CC  CO 

co  co  co 

CO  CO  CO 

30  33  CO 

coco  co 

CO  CO  CO 

OnOOn 

O'  O'  O' 

OOCT' 

O'  O'  O' 

rs  rr-  o 

O'  O'  O' 

'  4 

o 

■H 

-P 

-P 

CO  rH  S) 

rH  <AlA 

hhd 

S3  LA  CO 

S3  S3 -4 

r-om 

a3 

© 

co  r—  4 

CO  o  o 

OOffl 

4S3  O' 

O'  OOO 

S3  LA  cA 

> 

© 

co  co  co 

lAvO  S3 

S3  S3  la 

iAu\r- 

r-LA-4 

-44- -4 

© 

u 

<H 

v0^\O 

S3  S3  S3 

S3  S3  S3 

©©IS 

r-r— r- 

r-r-r- 

rH  <AIA 

cr  o  co 

_4  S3  r- 

O'CAcO 

cm  _4r- 

O' CM  LA 

© 

T3 

e 

CM  CO  CO 

CO  <A<A 

cArAcA 

A44 

LA  o  o 

O  rH  H 

3 

CM  CM  CM 

CM  CO  CM 

CO  CM  CM 

CM  CM  O 

Or — 1  rH 

rH  rH  rH 

-P 

rH  rH  rH 

rH  rH  rH 

i — 1  i — 1  a — i 

H  rl  H 

r — 1  rH  i — 1 

r — 1  rH  r — 1 

•H 

-P 

© 

CM  CM  CM 

CM  CM  CO 

CM  CM  CM 

CO  CM  CO 

CM  CM  CO 

CM  CM  CM 

o 

LALALA 

LALALA 

LA  LA  -A 

LALALA 

LALALA 

LALALA 

© 

T3 

_4  S3  r- 

O'  o  o 

rH  CO  CM 

CA-4G' 

O' 33  CO 

co  r—  r— 

-4-4-4 

fA-4-4 

-4-4-4 

-4.4  LA 

lALA  la 

LA  LA  LA 

3 

rH  rH  rH 

i — 1  < — 1  i — 1 

rH  rH  rH 

< — 1  i — 1  rH 

rH  rH  rH 

rH  rH  rH 

■P 

•H 

bO 

r-r-r- 

r-r-r- 

r—  r—  r— 

r-r-r- 

r—  r—  c— - 

r-r-r- 

C 

rH  rH  rH 

rH  rH  i — 1 

rH  rH  rH 

rH  rH  rH 

rH  t — 1  i — 1 

rH  rH  rH 

o 

o 

i — 1  i — 1  rH 

rH  rH  rH 

i — I  rH  rH 

i — 1  i — 1  r — 1 

rH  rH  rH 

rH  rH  rH 

c 

o 

CO  O'  o 

rH  CM  CA 

-4  LA  S3 

r-CO  rH 

CM  O  rH 

CM  ca_4 

-p 

4S 

CQ 

rH 

rH  rH 

rH  rH  rH 

Jx,  Pc 

O  C3 

CJ  O  CiJ 

C5  CS  PJ 

PI  PI  4 

PJ  PI  PI 

4 

o 

(A_4la 

S3  r—  CO 

O'  O  rH 

CM  CA-4 

LA  S3  r- 

CO  O'  o 

& 

r-r-r- 

r-r-r- 

r— co  co 

COCO  CO 

33  CO  CO 

coco  O' 
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P 

©  P 

© 

b-  -d-LA 

cOLA_d 

co  oao 

O  CM  O 

CO  b-LA 

Cb  Cb  i — 1 

G  © 

bO 

*  *  # 

4  4  4 

4  4  4 

4  4  4 

•  •  • 

4  4  4 

60  g 

tH 

r—  b-  a- 

CO  C — O 

LA  LA  CO 

cO  CM  P 

O  COLA 

co  CM  O 

G  o 

1 - 1 

Cb  Cb  Cb 

Cb  Cb  Cb 

Cb  Cb  Cb 

Cb  Cb  Cb 

Cbco  co 

CO  GO  CO 

o  c 

rH 

H  rl  H 

p  p  p 

P  P  P 

P  P  P 

p  p  p 

P  P  P 

CQ  < 

P 

G 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

4- 

OQ 

rH 

•H  j3  M 

© 

COOIA 

OLACb 

b-CO  CO 

O  P  o 

co_d  (b 

O  ajla 

©  o  l 

bn 

4  •  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

G  ©  < 

•H 

CO  O  CO 

CO  C —  O 

© 

b-  C —  b- 

'0-0  0 

b-b-b- 

U 

®  G 

Eh  o 

■g 

Li  - 

©  G  t>s 

m 

P 

© 

OOP 

COb-CO 

CO  LA  CO 

o  cb-4 

0-4  la 

O  cb  O 

P  ©  P 

bfl 

la -d  o 

rO  Cb  CO 

_d-co  O 

Cb  CM  O 

©Pd 

Cb  o  o 

cl  G  as 

«H 

4  •  • 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

6  MS 

rH 

O  'Ll  O 

© 

CM  P  O 

Cb  CbCO 

b-LA  AJ 

O  O  b- 

P  G  O 

P 

o  o  o 

o  o  o 

O  O  O 

Cb  Cb  Cb 

Cb  Cb  Cb 

cb  cbco 

CO  O  C 

P 

CM  CM  CM 

CM  CM  CM 

CM  CM  CM 

P  P  P 

Pi — IP 

Pi — IP 

CO  < 

G 

1  1  1 

i  i  i 

i  i  i 

1  1  » 

1  1  1 

1  1  1 

F. 

ro 

•H  t>> 

P 

fH 

© 

b-LA  CO 

-dO-d- 

-dCM  CM 

P  OLA 

Cbcob- 

c— P  o 

s 

bD 

O  rH  O 

b-_d-0 

P  Cb  P 

CO  CM  O 

CO  LA  O 

AJ  O  O 

©  G 

P 

•  ©  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

0  4  4 

©  o 
u  g 

rH 

fH 

©lAf^ 

-4-d*-4 

P  CM  co 

CO  CM  CO 
-d--4-4 

CO  co  CM 
-4-4  ~d- 

HOH 

P  O  Cb 
-4-4  co 

fe  <c 

•H 

G 

'Cf 

rH  H  O' 

r —  r — o 

CO  LA  LA 

COP  co 

LA  b-LA 

CO  AJ  CM 

© 

00 

LALA  AJ 

co  _4p 

_40  b- 

co_d-CM 

LA  Cb  LA 

LAO  b- 

>  -P 

rH 

4  •  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

© 

O  CM  f 

CO  P  LA 

b-  Cb  CM 
cOrO-d- 

4©  cb 
dr -d-d 

CO  CM  Cb 

cor- LA 

©  > 

bfl 

CM  CM  CM 

CM  coco 

LAO  O 

b-b- CO 

m  © 

O  g 

d 

O  O  O 

O  O  O 

O  O  O 

o  vO  o 

o  o  o 

o  o  o 

o  o 

d 

O  O  O 

o  o  o 

O  o  o 

o  6  o 

o  o  o 

O  O  o 

G 

co  coco 

co  coco 

30  CO  CO 

CO  to  co 

co  co  co 

COCO  CO 

C"  O'-  Cb 

Cb  Cb  Cb 

Cb  Cb  Cb 

Cb  cb  Cb 

Cb  Cb  Cb 

Cb  Cb  Cb 

~TT 

o 

•H 

•P 

p 

LALA-d* 

_4LA-d" 

CM  b-b- 

CO  OLA 

CM  O  b- 

b-_4co 

© 

© 

HCOrl 

CO  O  CO 

P  CO  LA 

-d  AJ  DO 

co  CM  LA 

p  b-b- 

> 

© 

4*  coco 

CM  CM  CM 

CM  P  P 

P  P  O 

O  CbCO 

CO  b-O 

© 

b-b-b- 

b-b-b- 

b-b-  b- 

b-C^b- 

b-O  o 

O  O  O 

rH 

CO  CM  Cb 

cOb-LA 

Cb  COO 

Cb O  CO 

b-b- CM 

LA  b-co 

© 

T5 

e 

P  CM  CM 

coco^t 

-d-LALA 

LAO  O 

O  P  CM 

CM  AJ  CO 

G 

p  p  p 

P  P  P 

p  p  p 

P  AJ  CM 

CM  AJ  CM 

AJ  CM  AJ 

4-> 

«=> 

p  p  p 

P  P  P 

p  p  p 

P  P  P 

P  P  P 

P  P  P 

P 

© 

CM  CM  CM 

CM  CM  CM 

CM  (M  CM 

AJ  AJ  AJ 

CM  AJ  CM 

AJ  AJ  AJ 

id 

o 

LA  LA  LA 

LA  LA  lA 

LA  LA  LA 

LA  LA  LA 

LA  LA  LA 

LAU\  -A 

N©  J- 

-4  cop 

O  O  Cb 

CbCO  CO 

b-LA  .4 

coco  AJ 

© 

•  •  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

4  4  4 

'd 

•"* 

LA  LALA 

LA  LA  LA 

LALA-d 

-j 

..rf  ,-f  -j: f 

G 

1) 

P  P  P 

p  p  1 — 1 

P  P  P 

p  p  p 

p  p  i — i 

• — i  p  p 

•H 

bO 

b~  r —  b— 

b-b-b- 

b-C^b- 

b-  b-b- 

b-b-b- 

b-b-b- 

G 

o 

p  p  p 

« — 1  1 — 1  1 — 1 

P  P  P 

P  P  P 

p  p  p 

p  p  p 

O 

p  p  p 

Pi — IP 

i — 1  P  P 

P  P  i — 1 

p  p  p 

p  p  p 

id 

G 

o 

•H 

LAOOO 

CbO  CM 

cO-d-LA 

o  b-co 

Cb  P  CO 

-dLAO 

p 

© 

p 

P  P  P 

P  CM  CM 

CM  CM  AJ 

AJ  AJ  CM 

CM  coco 

CO  co  co 

>d  4  id 

►d  id  .d 

.d  .d  4 

^q  »d  .d 

►d  »d  id 

•did  4 

CO 

• 

o 

PCMm 

-4  LAO 

b-CO  Cb 

O  P  CM 

co_4la 

O  b-co 

!3 

Cb  Cb  Cb 

Cb  Cb  Cb 

Cb  Cb  Cb 

O  O  O 

O  O  o 

o  o  o 

P  P  P 

p  p  p 

i — 1  P  i — 1 
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o 

< 

£ 

co 

a 

Eh 
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M 

O 

5 

o 

c 

M 

§ 

S 
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WJ 

rH 

©  rH 

© 

PLA_d 

A- LA  M3 

i — 1  LA  t — 1 

LA_d  CA 

OPP 

P  O 

P 

P  CD 

bo 

•  •  • 

•  •  ♦ 

m  •  • 

•  •  • 

•  •  • 

•  • 

9 

bO  £ 

■H 

LAIAO 

CA_c±_d 

A-P)  P 

P  rH  CA 

CAO  O 

P  P 

CA 

P  o 

1 - 1 

pp  ca 

CO  CA  CA 

CO  00  00 

GO  CO  CO 

CA  CAP 

OO  p 

P 

o  c 

CD  < 

rH 

H  H  H 

rH  rH  rH 

rH  rH  rH 

rH  rH  rH 

i — 1  r — 1  i — | 

rH‘rH 

1 — 1 

•H 

£ 

1  1  1 

1  1  1 

1  1  1 

1  1  1 

1  t  | 

i  i 

c§ 

•H  -H 
©  +?  H 
£  8  * 

co 

rH 

© 

bO 

•H 

p  p  cm 

O  LA  rH 

LAP)  P> 

PIA  CA 

-d-d- la 

OP 

P 

• 

H 

O' O' CM 

CM  CO  CO 

COCO  A- 

P  P_p" 

CM  rH  rH 

O  CA 

P 

©  L  ' — 

E— 1  5 

a 

rH 

1 - 1 

rH 

rH  rH  rH 

rH 

o 

£ 

co 

(D  U  t>a 

rH 

cd 

cAO  cA 

A-A-  O 

LA  CO  <A 

OP  CA 

P  o  CM 

A- CMP 

CA 

i — 1  ®  i — 1 

bO 

-otAJP 

P  CAA- 

P)  O  A- 

rH  CO  A- 

CAO  ca 

A-P  A- 

-d- 

Cl  P  © 

•H 

9 

£  bQ  £ 

rH 

LA  LA  CM 

rH  CM  CM 

LAVACA 

CAP  P 

LA  CM  rH 

P  LA  Al 

AJ 

•H  2  O 

rH 

COCO  O 

o  o  o 

CA  CA  CA 

CAP  CA 

O  O  O 

CA  CAP 

P 

woe 

•H 

H  H  CM 

CM  CM  CM 

rH  rH  rH 

rH  rH  rH 

CM  CM  CM 

Hr — 1  rH 

rH 

CD  <© 

£ 

1  1  1 

1  1  1 

I  1  1 

1  1  1 

1  1  1 

1  1  1 

i 

U 

w 

rH 

•H  f>5 

© 

rHP-d" 

P  CM  CA 

rH  O  CM 

CA  O  CA 

rH  P  P 

P  CA  O 

p 

<  rH 

bO 
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AHU\ 

CO  CM  A- 

p  CAP 

A-p_d- 
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p 

© 

•H 

•  •  • 

•  •  * 

•  •  • 

•  •  • 

•  •  • 

•  •  • 

• 

©  £ 

rH 

O'  o  o 
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CA  CM  CM 

LA  A-P 

P  CAA- 

rH 

©  o 
u  G 
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•H 
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la -d-d* 

LALA-d- 
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-d_d-^d 

-d-d-cA 
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£ 

-d^P  A- 

LA  LAP) 

P)  O  A- 

r — 1  CA  i — 1 

CM  CA_d 

CA  CA  CA 

'O 

© 

LA  O  CA 

P>  P±CM 

rH  A- LA 

p_d-o 

O  A-  CM 

CA  AJ  A- 

P 

©  r*s 

rH 

• 

>  -P 

© 

CACACA 

LAO  CM 

HLA  A- 

CM  A-  CM 

_d"CA_d 

AJ  LA  CA 

P 

bO 

O  O  cm 

CM  CAtA 

-d'-dLA 

P  A-  rH 

CM  CM  CA 

-dvr.  O 

CA 

©  > 

•H 

C —  P  p 

P  P  P 

p  p  p 

P  P  P 

P  P  P 

P  P  A- 

P 
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IRAIN  CORRECTIONS  FOR  HAYFORD  ZONES  B  TO  J  ON  THE  ATHABASKA  GLACIER 


Units  of  Milligals 


No. 

Station 

B 

c 

D 

E 

F 

G 

H 

I 

J 

Total 

1 

A 

1 

0.16 

0.19 

1.11 

3.57 

3.65 

1.76 

0.76 

0.35 

0.19 

11.74 

2 

A 

2 

0.20 

0.18 

0.83 

3.50 

3.70 

1.80 

0.80 

0.34 

0.18 

n.53 

3 

A 

3 

0.20 

0.19 

0.65 

3.44 

3.81 

1.90 

0.83 

0.33 

0.18 

11.53 

4 

A 

4 

0.20 

0.19 

0.37 

2.56 

3.65 

1.98 

0.87 

0.32 

0.18 

10.32 

5 

A 

5 

0.24; 

0.18 

0.30 

1.69 

3.51 

2.05 

0.92 

0.32 

0.18 

9.39 

6 

A 

6 

0.20 

0.16 

0.28 

1.55 

3.53 

2.10 

0.95 

0.32 

0*18 

9.27 

7 

A 

7 

0.21 

0.15 

0.30 

1.42 

3.56 

2.26 

0.97 

0.31 

0.18 

9.36 

8 

B 

1 

0.15 

0.26 

1.04 

3.06 

4.54 

2.32 

1.05 

0.49 

0.23 

13.14 

9 

B 

2 

0.00 

0.15 

0.95 

2.51 

4.41 

2.33 

1.07 

0.48 

0.23 

12.13 

10 

B 

3 

0.00 

0.08 

0.1)5 

1.96 

4.28 

2.34 

1.09 

0.47 

0.23 

10.90 

11 

B 

4 

0.00 

0.06 

0.22 

1.75 

4.09 

2.35 

1.14 

0.46 

0.23 

10.30 

12 

B 

5 

0.00 

0.06  0.14 

1.54 

3.93 

2.36 

1.18 

o.45 

0.23 

9.89 

13 

B 

6 

0.02 

0.06 

0.09 

1.21 

3.82 

2.37 

11.22  0.44 

0.23 

9.46 

14 

B 

7 

0.00 

0.06 

0.11 

1.32 

3.60 

2.38 

1.23 

0.43 

0.23 

9.36 

15 

B 

8 

0.00 

0.06 

0.26 

1.43 

3.39 

2.41 

1.25  0.42 

0.23 

9.45 

16 

B 

9 

0.03 

0.10 

0.80 

2.19 

2.80 

2.42 

lo32 

0.40 

0.23 

10.29 

17 

B 

10 

o.ei 

0.62 

1.60 

2.54 

2.43 

2.44 

1.36 

o.37 

0.23 

12.40 

18 

C 

1 

0.58 

0.84 

1.02 

2.85 

4.20 

2.43- 

1.02 

o.65 

0.27 

13.84 

19 

C 

2 

1.22 

10.40 

0.95 

2.46 

4.00 

2. ill 

1.08 

0.64 

0.27 

12.43 

20 

c 

3 

1.00 

0.04 

0.58 

2.08 

3.79 

2.42 

1.14 

0.63 

0.27 

!io.95 

21 

c 

4 

0.00 

0.06 

0.35 

1.74 

3.64 

2.43 

1.20 

0.62 

0.27 

10.31 

22 

c 

5 

0.00 

0.04 

0.33 

1.40 

3.49 

'2.43 

1.26 

10.61 

0.27 

9.83 

23 

c 

6 

0.00 

o.05 

0.05 

1.20 

3.34 

2.44 

1.32 

0.60 

0.27 

9.27 

24 

c 

7 

0.00 

1 

o.05 

0.02 

0.99 

3.19 

.2.44 

1.38 

0.59 

0.27 

8.93 

25 

c 

8 

0.00 

0.05 

0.02 

1.02 

2.77 

2.44  P..45 

0.58 

0.27 

8.60 

26 

c 

9 

O.oo 

0.04 

0.02 

1.04 

2.75 

I2.U4  1.52  0.57 

0.27 

!  8.65 

27 

c 

10 

0.00 

0.04 

jO. 06 

1.03 

2.67 

(2.45  fL .53 

0.56 

10.27 

8.61 

28 

c 

11 

0.00 

0.04 

0.09 

1.03 

2.63 

2.44  1.55 

b.55 

0.27 

8.60 

29 

c 

12 

0.00 

0.04 

0.29 

1.28 

2.17 

2.44  1.56 

0.54 

0.27 

8.59 

30 

c 

13 

0.00 

0.08 

o.44 

1.53 

2.14 

2.43 

1.58 

b  .53 

0.27 

9.00 

io  6 

APPENDIX  2  -  TABLE  III 

TERRAIN  CORRECTIONS  FOR  HAYFORD  ZONES  B  TO  J  ON  THE  ATHABASKA  GLACIER 

Units  of  Milligals 


No. 

Station 

B 

C 

D 

E 

F 

G 

H 

I 

J 

Total 

31 

C 

14 

0.02 

0.17 

1.30 

1.85 

1.95 

2.43 

1.59 

0.52 

0.27 

10.10 

32 

C 

14* 

0.03 

0.37 

1.87 

2.18 

1.76 

2.43 

1.61 

o.5o 

0.27 

11.02 

33 

C 

16 

0.56 

0.89 

2.54 

2.18 

1.76 

2.43 

1.61 

o.5o 

0.27 

12.74 

34 

D 

1 

0.33 

0.13 

0.87 

1.80 

3-08 

2.45 

0.98 

0.69 

0.32 

10.65 

35 

D 

2 

0.09 

0.27 

0.82 

1.70 

2.89 

2.44 

1.05 

0.70 

0.32 

10.28 

36 

D 

3 

0.00 

0.01 

0.43 

1.58 

2.76 

2.42 

1.12 

0.70 

0.32 

9.34 

37 

D 

4 

0.00 

0.03 

0.26 

1.21 

2.61 

2.40 

1.19 

0.70 

0.32 

8.72 

38 

D 

5 

0.00 

0.02 

0.08 

o.85 

2.45 

2.36 

1.26 

0.69 

0.32 

6.03 

39 

D 

6 

0.00 

0.02 

0.00 

0.85 

1.93 

2.30 

1.33 

0.69 

0.32 

7.44 

ko 

D 

7 

0.00 

0.02 

0.04 

0.86 

1-53 

2.25 

1.41 

0.69 

0.32 

7.12 

41 

D 

8 

0.00 

0.02 

0.07 

0.94 

1.54 

2.24 

1.44 

0 . 68 

0.32 

7.25 

42 

D 

9 

0.00 

0.02 

0.16 

1.03 

1.58 

2.22 

1.47 

0.67 

0.32 

7.47 

43 

D 

10 

0.00 

0.01 

0.35 

1.06 

1.47 

2.20 

1.51 

0.66 

0.32 

7.58 

44 

D 

11 

0.00 

0.10 

1.02 

1.09 

1.36 

2.19 

1.54 

0.66 

0.32 

8.28 

45 

D 

lit 

0.03 

0.63 

1.78 

0.84 

1.04 

2.17 

1.58 

0.65 

0.32 

9.04 

46 

E 

0 

0.23 

0.16 

0.86 

1.96 

2.45 

2.53 

1.29 

0.72 

0.37 

10.57 

47 

E 

1 

0.10 

0.27 

0.75 

2.00 

2.32 

2.45 

1.29 

0.75 

0.37 

10.30 

48 

E 

2 

0.13 

0.09 

0.57 

2.04 

2.19 

2.36 

1.29 

0.78 

0.37 

9.82 

49 

E 

3 

0.08 

0.10 

0.26 

1.33 

2.13 

2.27 

1.29 

0.81 

0.37 

8.64 

50 

E 

4 

0.08 

0.08 

0.07 

0.62 

2.08 

2.18 

1.30 

0.84 

0.37 

7.62 

51 

E 

5 

0.12 

0.08 

0.09 

0.50 

1.86 

2.0  Q 

1.30 

0.87 

0.37 

7.28 

52 

E 

6 

0.14 

0.09 

0.10 

0.68 

1.39 

2.03 

1.38 

0.88 

0.38 

7.07 

53 

E 

7 

0.10 

0.14 

0.09 

0.67 

1.26 

1.96 

1.40 

0.90 

0.38 

6.90 

54 

E 

8 

0.00 

0.13 

0.28 

0.66 

1.13 

1.89 

1.50 

0.93 

0.38 

6.90 

55 

E 

9-C 

0.07 

0.43 

0.57 

0.71 

o.84 

1.82 

1.62 

0.96 

0.38 

7.40 

56 

H 

1 

0.29 

0.19 

0.51 

1.21 

2.21 

2.33 

1.27 

0.96 

0.42 

9.39 

57 

H 

2 

0.15 

0.18 

0.41 

0.95 

2.00 

2.28 

1.32 

0.97 

0.42 

8.68 

58 

H 

3 

0.06 

0.10 

0.10 

0.53 

1.61 

2.14 

1.42 

1.01 

0.42 

7.39 

59 

H 

4 

0.00 

0.08 

0.07 

0.35 

1.33 

1.94 

1.52 

1.02 

0.42 

6.73 

60 

H 

5 

0.00 

0.10 

0.10 

0.37 

1.16 

1.82 

1.55 

1.03 

0.42 

6.55 

' 
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TERRAIN  CORRECTIONS  FOR  HAYFORD  ZONES  B  TO  J  ON  THE  ATHABASKA  GLACIER 

Units  of  Milligals 


No. 

Statior 

B 

C 

D 

E 

F 

G 

H 

I 

J 

Total 

61 

H 

6 

0.00 

0.10 

0.15 

0.39 

0.99 

1.76 

1.59 

1.03 

0.42 

6.43 

62 

H 

7 

0.00 

0.07 

0.17 

0.40 

0.90 

1.73 

1.59 

1.03 

0.42 

6*31 

63 

H 

8 

0.00 

0.19 

0.29 

0.43 

0.80 

1.68 

1.60 

1.03 

0.42 

6*44 

64 

H 

9 

0.1+3 

0.25 

0.60 

0.46 

0.60 

1.61 

1.60 

1.03 

0.42 

7.00 

65 

F 

0 

0.35 

0.22 

0.41 

0.90, 

1.95 

2.14 

1.33 

1.03 

0.46 

6.79 

66 

F 

1 

0.53 

0.19 

0.29 

0.78 

1.77 

2*02 

1.36 

1.15 

0.46 

6.55 

67 

F 

2 

0.27 

0.10 

0.18 

0.63 

1.64 

2.00 

1.40 

1.16 

0.46 

7.84 

66 

F 

3 

0.14 

0.09 

0.15 

0.48 

1.51 

1.98 

1.45 

1.16 

0.46 

7.42 

69 

F 

4 

0.00 

0.14 

0.13 

0.37 

1.35 

1.92 

1.52 

1.16 

0*46 

7.05 

70 

F 

5 

0.00 

0.12 

0.17 

0.27 

1.18 

1.81 

1.58 

1.16 

0.46 

fc.75 

71 

F 

6 

0.00 

0.15 

0.23 

0.24 

1.08 

1.75 

1.60 

1.17 

0.46 

6.68 

72 

F 

7 

0.00 

0.15 

0.25 

0.22 

1.02 

1.66 

1.63 

1.19 

0.46 

6.58 

73 

F 

8 

0.00 

0.15 

0.25 

0.22 

1.02 

1.66 

1.63 

1.19 

0.46 

6.58 

74 

F 

9 

0.24 

0.17 

0.29 

0.40 

0.75 

1.63 

1.74 

1.22 

0.45 

6.89 

75 

F 

10 

0.34 

0.30 

0.35 

o.44 

0.70 

1.64 

1.79 

1.24 

0.45 

7.25 

76 

G 

1 

0.22 

0.21 

0.22 

0.76 

1.56 

1.83 

2.05 

1.37 

0.52 

8.74 

77 

G 

2 

p.25 

0.21 

0.30 

0.58 

1.41 

1.80 

2.05 

1.37 

0*52 

8.49 

78 

G 

3 

0.23 

0.24 

0.36 

0.41 

1.26 

1.75 

2.05 

1.37 

0*52 

8.19 

79 

G 

4 

0.17 

0.24 

0.35 

0.35 

1.16 

1.69 

2.06 

1.36 

0.52 

7.90 

80 

G 

5 

p.31 

0.36 

0.J9 

0,30 

1.06 

1.64 

2.06 

1.35 

0.52 

7.99 

81 

G 

6 

0.26 

0.36 

0.39 

0.24 

1.02 

1.59 

2.07 

1.35 

0.52 

7.80 

82 

G 

7 

0.40 

0.46 

0.37 

0.22 

1.00 

1.59 

2.06 

1.34 

0.62 

7.96 

83 

G 

8 

0.46 

0.78 

0.59 

0.19 

0.99 

1.58 

2.06 

1.34 

0.62 

8*61 

84 

L 

1 

p.o4 

0.05 

0.47 

2.64 

3.87 

1.95 

0.85 

0.30 

0*17 

10.34 

85 

L 

2 

b.io 

0.12 

0.27 

2.10 

3.70 

1.99 

0.88 

0.32 

0.17 

6.65 

86 

L 

10 

0.18 

0.09 

0.17 

1.40 

4.08 

2.25 

1.06 

0.38 

0.2i 

9.82 

87 

L 

11 

p.05 

0.07 

0.16 

1.30 

3.90 

2.30 

1.13 

0*39 

0.21 

9.51 

88 

L 

12 

0.08 

0.04 

0.16 

1.21 

3.82 

2.38 

1.22 

0.42 

0.22 

Q.55 

89 

L 

13 

0.00 

0.04 

0.13 

1.13 

3.60 

2.40 

1.25 

0.44 

0.23 

0 .22 

90 

L 

14 

0.00 

! 

0.04 

0.07 

1.07 

3.40 

2.41 

1.27 

0.41 

0.24 

6.91 
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TERRAIN  CORRECTIONS  FOR  HAYFORD  ZONES  B  TO  J  ON  THE  ATHABASKA  GLACIER 

Units  of  Milligals 


No. 

Station 

B 

c 

D 

E 

p 

G 

H 

I 

J 

Total 

91 

L  15 

0.00 

0.06 

0.06 

1.01 

3.19 

2.42 

1.32 

0.48 

0.25 

8.79 

92 

L  16 

0.00 

0.11 

0.03 

1.02 

3.17 

2.43 

1.41 

0.53 

0.26 

8.96 

93 

L  18 

0.00 

0.04 

0.01 

0.95 

2.70 

2.43 

1.49 

0.60 

0.28 

8.50 

94 

L  19 

0.00 

0.03 

0.01 

0.86 

2.27 

2.42 

1.46 

0.63 

0.29 

7.97 

93 

L  20 

0.00 

0.02 

0.01 

0.86 

0.88 

2.35 

1.42 

0.66 

0.31 

7.51 

96 

L  22 

0.00 

0.02 

0.01 

0.73 

1.46 

2.23 

1.37 

0.72 

0.33 

6.87 

97 

L  23 

0.00 

0.03 

0.04 

0.63 

1.40 

2.21 

1.32 

0.76 

0.35 

6.74 

98 

L  24 

0.00 

0.05 

0.05 

0.65 

i.4o 

2.20 

1.30 

0.78 

0.36 

6.79 

99 

L  25 

0.00 

0.06 

0.13 

0.62 

1.40 

2.11 

1.30 

0.83 

0.37 

6.82 

100 

L  26 

0.10 

0.08 

0.10 

0.68 

i.4o 

2.09 

1.30 

0.87 

0.37 

6.99 

101 

L  27 

0.16 

0.14 

0.08 

0.68 

1.39 

2.00 

1.38 

0.88 

0.38 

7.09 

102 

L  28 

0.14 

0.11 

0.13 

0.60 

1.30 

1.97 

1.39 

0.95 

0.39 

6.98 

103 

L  29 

0.14 

0.12 

0.12 

c.50 

1.15 

1.90 

i.44 
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6.76 

104 
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1.00 
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1.60 
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o.44 

6.44 

io5 

L  33 
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1.00 

1.68 

1.68 

1.23 

o.47 

6.89 

106 

L  34 
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0.26 

1.00 

1.60 

1.75 

1.26 
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7*03 

107 

L  35 

p.18 

0.21 
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0.30 
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1.58 

1.81 

1.27 

0.49 

7.23 

108 

L  36 
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0.26 

1.00 

1.58 

1.98 

1.33 

o.5i 

7.51 

109 

Lake 

0.00 
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1.66 

1.93 

2.54 

1.68 

0.77 

9.73 

110 

Bridge 

0.00 

0.06 

0.24 

o.84 

1.66 

1 .63 

2.54 

1.68 

0.77 

9.72 

111 

M  13 

0.36 

0.69 

1.73 

2.60 

2.30 

2.44 

1.42 

0.41 

0.23 

12.18 

112 

M  14 

0.42 

0.78 

1.95 

2.10 

2.20 

2.44 

1.48 

0.43 

0.24 

12.04 

113 

M  18 

0.00 

0.03 

0.49 

1.36 

1.80 

2.40 

1.55 

o.57 

0.28 

8.48 

114 

M  19 

0.00 

0.04 

0.54 

1.21 

1.49 

2.35 

1.55 

0.61 

0.29 

8.08 

115 

M  23 

o..  23 

0.49 

1.45 

0.75 

0.90 

2.10 

1.55 

0.73 

0.34 

8.54 

116 

M  24 

p.48 

0.47 

1.52 

0.67 

0.79 

2.00 

1.57 

0.79 

0.36 

8.65 

117 

M  28 

0 . 65 

0.45 

0.38 

0.68 

0.70 

1.74 

1.60 

1.01 

0.41 

7.62 

118 

M  31 

0.22 

0.18 

0.35 

0.45 

0.70 

1.63 

1.62 

1.05 

0.43 

6.63 

119 

M  35 

0.59 

0.40 

0.46 

0.44 

0.49 

1.57 

1.85 

1.24 

0.47 

7.51 

120 

Gravel 

1.31 

0.34 

1.48 

4.26 

3.84 

2.10 

0.92 

0.45 

0.20 

14. 9< 
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TERRAIN  CORRECTIONS  FOR  HAYFORD  ZONES  B  TO  J  ON  T!JE  ATFABASKA  GLACIER 

Units  of  Milligals 


No. 

Station 

B 

C 

D 

E 

F 

G 

H 

I 

J 

Total 

121 

B 

0.07 

0.22 

1.15 

2.48 

4.19 

2.39 

1.05 

0.56 

0.26 

12.37 

122 

D 

0.46 

0.42 

0.76 

1.96 

3.47 

2.14-2 

1.01 

0.66 

0.28 

11.44 

123 

E 

0.62 

o.Jt-5 

0.96 

1.88 

3.20 

2.44 

0.98 

0.68 

0.31 

11.52 

124 

G 

0.25 

0.29 

0.72 

1.60 

2.41 

2.47 

1.15 

0.76 

0.34 

9.99 

125 

H 

0.27 

0.13 

0.55 

1.38 

2,29 

2.40 

1.29 

0.92 

0.39 

9.62 

fT 

APPENDIX  3  -  TABLE  VI 

RESULTS  OP  COMPUTER  PROGRAM  FOR  THE  ATFABASKa  GLACIER 
Gravity  anomaly  in  units  of  milligals:  Depth  in  feet. 


Line 

6666 

Station 

Model  I 

Model  II 

Model  III 

Maximum  Depth 

I,  II,  III 

A 

A 

A 

A 

1 

3 

5 

7 

I}-. 90 
9.14 
12.03 
11.60 

3.91+ 

7.21+ 

lO.ltO 

10.35 

700,  560  -- 

B 

l 

10.82 

9.82 

9.93 

1,070,  860,  1( 

B 

3 

13.61 

12.33 

13.3lt 

L 

12 

15.13 

13.51 

15.34 

B 

9 

10.87 

9.78 

12.23 

B 

10 

5.H 

b.5b 

7.19 

C 

0 

3.35 

2.65 

1020,  820 

c 

3 

11.62 

9.99 

c 

5 

lit.  35 

12.54 

c 

7 

15.76 

13.71 

c 

9 

15.86 

13.73 

c 

12 

13-1+3 

11.51 

c 

lb 

9.72 

8.29 

c 

16 

3.68 

2.99 

D 

1 

8.97 

7. 47 

8.20 

1000,  800,  10( 

D 

2 

11.32 

9.28 

10.51 

D 

k 

lit.  30 

11.93 

13.90 

I 

D 

6 

15.1+3 

13.20 

15.37 

D 

10, 

12.93 

10.72 

12.94 

D 

hi 

6.92 

5.81 

7.65 

E 

0 

5.88 

5.  ok 

800,  640 

E 

2 

8.88 

7.83 

E 

3 

11.82 

10.14 

E 

5 

1  13.37 

n.64 

E 

7 

12.46 

10.70 

E 

9C 

- - 

7.23 

J - 

6.16 

- - 

*  ■  3K  j.  3ffH  "t,: *  I  1 
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RESULTS  OF  COMPUTER  PROGRAM  FOR  E  ATHARASKA  GLACIER 


Gravity  anomaly  in  units  of  milligals:  Depth  in  feet. 


Line 

6666 

Station 

. 

Model  I 

Model  II 

Model  III 

Maximum  Depth 

I,  11^  III 

H 

H 

H 

H 

H 
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11.56 
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6.68 
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7.63 
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8.69 
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6.63 
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10 

14-.13 

3.22 

.  .  .  
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2.48 
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3.78 

6.06 

2.70 
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6 
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5.73 
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1 
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FIGURE  3 

BOUGUER  GRAVITY  MAP 

of  the 

ATHABASKA  GLACIER 

ALBERTA  CANADA 
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